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lattice theory a recent prediction for the (logarithmic) decay of lattice artifacts 
is probed. 
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1. Introduction 



The XY-model in two dimensions is of prime interest in the field of statistical mechanics, 
intriguing in particular because of its unusual phase transition. On general grounds one 
expects that a suitable scaling limit in the high temperature phase gives rise to a massive 
relativistic quantum field theory (QFT). Though an enormous literature exists on the 
statistical mechanics aspects of the system, to the best of our knowledge the nature of 
this QFT has never been systematically explored. As part of a long term project on 
quantum non-linear sigma models we thus address here the question: 

What are the qualitative and quantitative features of the QFT obtained from the XY-model 
by taking the massive continuum limit? 

The question is of interest, first in that it highlights the important problem of controlling 
the approach to the continuum in the lattice formulation of QFTs, and second because 
the proposed continuum QFT seems to possess a rich, hitherto unknown superselection 
structure. The problem of controlling the continuum limit of a lattice system based on 
numerical simulations alone is notoriously difficult for several reasons, (i) First because 
one often lacks rigorous knowledge of the phase structure and the position(s) of the critical 
points where the correlation length becomes infinite, (ii) Secondly there are in general 
no rigorous results on how a quantity approaches its continuum limit as a function of the 
correlation length, even if the existence of the limit is taken for granted. 

Knowledge of (i) and (ii) is crucial not only for quantitative aspects but also for matters 
of principle, like the status of asymptotic freedom beyond perturbation theory in non- 
abelian models, or more generally the physical differences or similarities of the continuum 
limits of (spin) systems with abelian and non-abelian symmetries. So far the application 
of lattice techniques to the extraction of physical quantities has concentrated on the non- 
abelian models. This suffers unfortunately from both of the before-mentioned problems 
(i) and (ii). Concerning (ii), the usual working hypothesis, attributed to Symanzik, is 
that in (perturbatively) asymptotically free theories, physical quantities approach their 
continuum limit rather rapidly with power- like corrections 1/^^, with p a positive integer 
(up to multiplicative logarithmic corrections). 

On the other hand for the XY-model we do have rigorous information on point (i). In 
particular the model with standard action is known to have two phases, one massless and 
the other massive [17]. The order of the phase transition has been argued by Kosterlitz 
and Thouless (KT) [28] to be infinite and this picture has been supported by various 
numerical studies [19]. Moreover concerning point (ii), one of the present authors (J.B.) 
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[4] , has argued that for a certain class of lattice actions and for certain observables (like 
the S-matrix or a current two-point function) leading lattice artifacts do not depend on 
the choice of lattice action and are calculable. However, they vanish extremely slowly, 
generically as inverse powers of the logarithm of the correlation length ^ (e.g. ~ 1/ln^^). 
One of the goals of the present paper is to test this proposal through extensive numerical 
simulations. Its derivation takes advantage of the Sine-Gordon description of the KT 
transition introduced by Amit et. al. [1]. Applying a series of (not entirely rigorous) steps 
invoking universality and making a sequence of mappings their analysis also entails the 
tentative identification: 

The massive continuum limit of the XY model is related to the Sine- Gordon (SG) model 
at its extremal fixed point jdsc = V^, in that both systems share subsets of fields with 
identical correlation functions. 

Examples of shared fields are the Noether current (proportional to the dual of the gradient 
of the SG scalar) and the energy momentum tensor. Similar to Coleman's SG-Thirring 
correspondence [13, 26] the mapping between the fields does not preserve locality in 
general and is likely not to be strictly one-to-one. Nevertheless the correspondence is very 
useful because the SG model is integrable and for such systems a direct non-perturbative 
continuum approach exists to construct the QFT, referred to as the form factor bootstrap. 
In fact, the SG model is the prototype integrable model and it has played an important 
role in the development of the form factor bootstrap method. Its bootstrap S-matrix was 
proposed in [52] and a large part of Smirnov's book [43] is devoted to the study of its 
soliton form factors, where also the form factors of the SG scalar and the Noether current 
are given. The second purpose of the present paper is to initiate a bootstrap construction 
of the 0(2) model along similar lines. Since the form factor approach is largely blind with 
respect to the local structure we can borrow many mathematical techniques from the SG 
model, but the interpretation as form factors of certain local 0(2) quantum fields requires 
careful justification, one strategy being the comparison with lattice simulations. To have 
a handy terminology we shall refer to the QFT defined through the massive continuum 
limit of the lattice XY model as the "XY QFT" , and to the QFT defined via the form 
factor bootstrap as the "bootstrap 0(2) model" . The basic proposition to be tested is 
that both QFTs coincide. 

Generally the problem of operator identification and classification in the bootstrap frame- 
work rests on conserved quantum numbers. By definition the 0(2) model has a manifest 
0(2) symmetry. Remarkably on the level of the bootstrap S-matrix and the scattering 
states a symmetry enhancement takes place: They are covariant with respect to a larger 
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non-abelian quantum group symmetry. It turns out that depending on the nature of cer- 
tain statistics phases the functional equations characterizing the form factors only have 
the manifest 0(2) symmetry or are covariant with respect to the hidden quantum group 
symmetry. This implies that the field operators of the bootstrap 0(2) QFT fall into two 
classes: those that are (trivial) 0(2) multiplets and those that are members of a nontrivial 
quantum group multiplct. For example a complete set of scattering states seems to be 
generated both by the spin field and a /ocaZ parafermion field of Lorentz spin 1/4. Both are 
relatively nonlocal and the latter is quantum group covariant while the former is not. Also 
the (one-component) Noether current of the 0(2) model is a member of a hidden isospin 
1 quantum group triplet, where however the charge ±2 partners can (already classically) 
not be expressed as local functions of the spin field. The energy-momentum tensor is a 
quantum group singlet. In view of the parafermion the 0(2) model possesses at least two 
super selection sectors; the full super selection structure and its relation to the quantum 
group multiplets remains to be explored. 

The lay-out of the paper is as follows. We start by briefiy recording the quantities con- 
sidered and introduce the bootstrap and the lattice formulation. In Section 3 we describe 
in more detail the bootstrap 0(2) model and discuss its quantum group invariance. We 
proceed with formulating the form factor equations and determine the statistics phases 
for which they exhibit the enhanced quantum group symmetry. Next the n < A particle 
candidate form factors of the spin field and the Noether current are obtained. Some de- 
tails on the quantum group structure and the form factor computation are relegated to 
appendices. 

The subsequent sections all involve lattice simulations, concerning which some general 
information is first collected in Section 4. Next, we report on measurements of the XY 
phase shifts using a refinement of the finite size technique first employed by Liischer and 
Wolff [31] for the 0(3) model. In Section 6 our measurements of the two-point functions of 
the spin field and the Noether current, and of the renormalized zero momentum coupling 
giR are presented. There we also compare the results with those obtained by the form factor 
computations, improving on our earlier estimate [7] . Finally in Section 7 we attempt some 
conclusions. 
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2. Bootstrap and Lattice 0(2) model 



One way to approach the 0(2) model is as the n = 2 member of the family of 0(n) 
nonlinear cr- mo dels with the classical Lagrangean 

C'^^'^^^d^S^d^S^ 5"5'^ = 1, a = l,2,...n. (2.1) 

From the viewpoint of classical field theory it may be surprising that this Lagrangean 
should correspond to a nontrivial QFT, since by substituting = cos0, 5*2 = sin0 it 
becomes quadratic, corresponding to a free theory. Also perturbatively the beta-function 
of the coupling g"^ vanishes identically for n = 2, suggesting a trivial scale invariant theory. 
The situation changes, however, when one is trying to nonperturbatively construct a QFT 
corresponding to (2.1). In this paper two such approaches are studied, both of which 
lead to a nontrivial massive QFT: the lattice approach which allows the construction 
of a massive continuum limit and the form factor bootstrap construction based on the 
indicated relation to the Sine- Gordon theory. These two constructions as well as the 
comparison of the resulting theories are the main subject of this investigation. 



2.1 Quantities to be investigated 

Clearly only physical quantities should be considered in this comparison. Wc shall study 
the S-matrix, the spin and current two-point functions, and the intrinsic coupling. The 
S-matrix will be discussed in sections 2.2 and 5; for the other quantities we collect here the 
key definitions. They apply to both formulations, and in fact to any 0(2) invariant scalar 
relativistic QFT with a mass gap. Let thus S'^{x) ,a = 1,2, denote a two-component 
(renormalized) scalar field (the "spin field"). For the Fourier transform of its Euclidean 
two-point function we write 

G{k)S''"'' = j (fxe'''''{S'''{x)S'''{0)) . (2.2) 

Its inverse is supposed to have the usual small momentum expansion 

G{k)-'^Z^'(^Ml + k' + 0{k')). (2.3) 

The coefficients can be expressed in terms of moments of the spectral density p(/i) via 

r, _ yll Ml _ 72 
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where M is the mass gap and 72 and 62 are the moments 



j^^M'Jdfi^, S^^M'Jdfi^. (2.5) 
Our normahzation for the spectral density is such that 

with the one-particle contribution given by p'-^^(/i) = S{n — M). (To avoid irrelevant 
complications we assume that the spectrum of the theory contains a doublet of stable 
particles of mass M.) 

The intrinsic or rcnormalized 4-point coupling is an important measure for the interaction 
strength of a QFT. A conventional definition is 



where G"'^'^'^ is defined through the Fourier transform of the connected 4-point function: 

4 

= {2Ti f5^^\ki + k2 + h + A;4)G'«i«2«3«^(A;i, k2, fcg, k^) . (2.8) 
The coupling g^i can then be written as 

»H = (2.9) 

where 74 is defined through 

^01020304^0 g g _ ^ ^a^a^ _|_ ^0103^0204 _|_ ^0104^0203 ^2 \Q'j 

In ref. [7] we computed the moments and coupling within the form factor approach in a 
certain truncation, and in Appendix D we present an improved approximation. 
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We also consider the two-point function of the Noether current J^: 



J d'xe''^J,ix)MQ)) = C6,, + ^ {k,K - k%.) • (2.11) 

Here C is the (regularization dependent) coefficient of a possible contact term. Only 
the coefficient I{k) of the transversal part is physical. It vanishes at zero momentum 
due to the assumed mass gap. The infinite momentum limit on the other hand is model 
dependent and can be finite or infinite. For the 0(2) model it is determined in Section 2.2. 



2.2 Integrability and bootstrap S-matrix 

A first hint why the 0(2) model (with classical Lagrangian (2.1) for n = 2) might be 
quantum integrable stems from the observation that the known bootstrap S-matrix of the 
0(n), n > 3, models has a smooth n ^ 2 limit [51]. Here we record this limit and also 
outline the relation to other integrable models. 

Assuming the spectrum of the model consisted of an 0(n) vector multiplet of massive 
particles the exact S-matrix of the n > 3 models was found by bootstrap methods [52]. 
For later use we adopt the projector decomposition 



\cd 



cd 



\cd 



(2.12) 



where 




_ {n - 2)0 + 27rt 
^'^^^ - {n- 2)9 - 27r. ^'^^^ ' 

-2tt 

, c ' -t- e 

■ sm (jjt 



1 + e- 



(2.13) 



The projectors are those on the 0(n) singlet, vector, and symmetric traceless tensors, i.e. 



cd 
ab 



_^ab^cd ^ 

n 

- ( 5'"'5^'^ 



(2.14) 



Contact to the Lagrangian (2.1) can be made through quantum conserved charges of 
higher spin that prevent particle production. Under mild extra assumptions Polyakov [41] 
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and Liischer [33] have shown the existence of such respectively local and nonlocal higher 
spin conserved charges. The latter in particular anticipate [33] a Yangian structure and 
entail the factorization equations that dictate the S-matrix. 

Much less is known about the 0(2) model. A simple observation is that the amplitudes 
(2.13) have a smooth n — 2 limit. This suggests that the 0(2) model might likewise 
be integrable and that its spectrum consists of a single 0(2) doublet of massive particles 
whose scattering is described by the n ^ 2 limit of the S-matrix (2.12). Although taking 
this formal n ^ 2 limit is not convincing in itself, the conclusion is corroborated by the 
KT theory [28] of the XY model and its reformulation in the context of the Sine-Gordon 
theory [1]. Before turning to the KT theory we thus briefly digress on the Sine-Gordon 
(SG) model. Its Lagrangean can be written as 

£SG ^ lo^^o^^ + ^ ~ ' ^2.15) 

where a has mass dimension 2 and /3 is the dimensionless SG coupling. It is also integrable 
and its spectrum and S-matrix was also found in [52]. The spectrum depends on (3 in 
a complicated way but it becomes simple for the range Stt > > An when it is free 
of bound states and consists of a single 0(2) vector of massive particles. In terms of 
z/ = || — 1 this corresponds to < < 1, and the S-matrix in this range can be written 
in the projector form (2.12), with n — 2, where now 

, , , sh^(i7r + e) , , , , , , ch^(i7r + ^) , , , 

Note that in the — > Stt (i^ — > 0) limit the SG S-matrix coincides with the n — > 2 limit 
of the 0{ri) S-matrix. 

Finally, in the vicinity of = Stt the SG model can also be related to a fermionic 
model [8] formulated in terms of a two-component Dirac fermion. It has a manifest SU(2) 
symmetry and is a variant of the chiral Gross-Neveu model with four-fermion interaction. 
The existence of this fermionic model sheds some light on the symmetry enhancement in 
the 0(2) model discussed in Section 3. For the details, however, the difference between 
SU-i{2) and SU{2) is crucial; c.f. the discussion at the end of Section 3.2. 

The relation to these other integrable QFTs can in particular be used to determine the in- 
finite momentum limit of the current two-point function I[k) in (2.11). In the 0(2) model 
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the limit is finite and exactly calculable. One way of computing /(oo) is by noting that it 
coincides with the coefficient of the Schwinger term in the current-current commutator. 
This commutator can be evaluated in the SG language using canonical quantization, and 
yields 

/(oo) = -. (2.17) 

TT 

The same result can be obtained using the relation to the before mentioned two-fermion 
model. Here, referring to the asymptotic freedom of the model in the fermion coupling 
constant, only a simple free fermion calculation has to be done. 



2.3 Standaird lattice action and KT theory 

The standard lattice action of the XY model is 

Sxy^kJ2[^- cos {^(x) - ip{x + jl))] . (2.18) 

X,IJ. 

We denote the inverse temperature (inverse of the bare couphng) of the XY model by K 
to avoid confusion with the SG coupling /3. 

This model has a high temperature phase at small K with exponential decay of correla- 
tions; it has been shown rigorously [17] that at low temperature (large K) the correlations 
decay only like a power of the distance. The exponential decay disappears therefore at 
a finite critical value Kc, this is the famous KT transition predicted by Kosterlitz and 
Thouless [28]. They argued that at not too small values of K typical configurations of 
the model can be described as a combination of 'smooth', topologically trivial configu- 
rations (spin waves) and a gas of vortices (of integer topological charge). The vortices 
have a logarithmic interaction and therefore form essentially a two-dimensional Coulomb 
gas, which has a transition from a high temperature phase with Debye screening to a low 
temperature dipole phase without screening. In the KT picture the transition is therefore 
described as 'vortex condensation'. A different perspective of this kind of phase transition 
was proposed in [38] according to which it is driven by the change from instanton-like 
defects (vortices) to super-instantons dominating at low temperatures. 

Prohhch and Spencer [17] established a rigorous version of the correspondence between 
the XY model and a type of Coulomb gas and used rigorous arguments inspired by the 
renormalization group (grouping of charges into neutral 'molecules') to show the absence 
of screening in this gas at low temperature. 
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Kosterlitz and Thouless employed heuristic energy-entropy considerations to show that 
in the transition region only vortices of topological charge ±1 are important and higher 
vortices can be neglected. It is easy to see that this system (spin waves and Coulomb gas 
with unit charge vortices only) is exactly equivalent to the SG model. In ref. [1] it was 
argued that the extremal SG fixed point (5* = -s/Stt , a* = is appropriate to describe 
the KT phase transition. The renormahzability of the SG model around this point was 
explicitly demonstrated up to two- loop order in a double expansion in a and 5 — ^ 



3. Bootstrap description and symmetry enhancement 

Here we detail on the proposed bootstrap S-matrix, the associated quantum group struc- 
ture and its implications for the form factors and the operator classification. 



3.1 S'C7_i(2) invciriance of the S-matrix 

The candidate S-matrix for the XY QFT can be rewritten as 



e 



\2 2m) ^ \2m) 



The S-matrix (3.1) satisfies the usual S-matrix postulates with the charge conjugation 
matrix Cab — ^ab and normalization S'^(O) = —^t^l- Note the nontrivial hmit 

5„t(±oo) = 5% - SabS^ , 5(±oo)^ = 1 . (3.2) 

The symmetry group of the massive 0(2) model is of course 0(2), as far as the La- 
grangian and the functional measure is concerned. The proposal (3.1) however entails 
that on the level of the S-matrix and the scattering states a symmetry enhancement 
takes place, in that on them a nonabelian symmetry operates. In view of the known 
Uq{su{2)) quantum group symmetry of the Sine-Gordon S-matrix [42] and the identifica- 
tion q = — e*'^(~^+^'^/^^\ one expects the symmetry to be W_i(s-u(2)). As a Lie algebra 
this is the same as Wi(sm(2)) = su{2), but the comultiplication in W_i(s'u(2)) differs from 
that in su{2). We refer to Appendix A for some basic definitions and our conventions on 
Uq{su{2)). To simplify the notation we shall write 5't/_i(2) for U-i{su{2)) from now on. 
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The easiest way to see the SU-i{2) invariance of the S-matrix (3.1) is to perform a 
projector decomposition. Defining S'^{6) := S^^{9) (so that S{0) = —1) it takes the form 

'iTT + 9 



s{e) = S2{e) 
1 



m-9 



Po + Pi 



with 



cd 
ab 



(3.3) 



Here PqPi — — PiPq and Pq + Pi — TL. Moreover Pq and Pi are the projectors onto the 
irreducible singlet and triplet representation of 5'C/_i(2), respectively. For comparison let 
us note that the SU{2) invariant S-matrix can likewise be written in the form (3.3), but 
the projectors are given by 

SU{2) : {P,r^ = i {S:St - 6^6^) , {Pi)li = \ {5:5t + 5t5t) . (3.4) 

The S-matrices (3.1) and (3.4) are of course also invariant under the usual real 0(2) 
transformations. It is often advantageous to diagonalize this action by means of a unitary 
basis transformation 




U 



(3.5) 



where we use Greek letters a, /3, . . . G {±} to label the components in the new basis. Ex- 
plicitly S]^ij{9) := U^U^S2f{9)U^U^, which now has Cap = Sa+i3fi as its charge conjugation 
matrix. Written in matrix form one finds the familiar pattern for (3.3) 



S{9) 



/ S2 \ 

St Sr 

Sr St 

V ^2 / 



St{9) 



S2{9), Sr{9)^-!^S2{9), (3.6) 



m — 9 



m — 9 



where the rows and columns refer to the (+-|-, H — , — h, ) ordering. For the SU{2) 

invariant S-matrix only the sign of St would be flipped, which in view of the previous 
discussion however indicates a very different group theoretical structure. Concretely the 
SU_i{2) invariance of (3.6) amounts to 



E± S{9) = S{9) E± , with Eh 










\ 




-1 









1 






v 





-1 1 





(3.7) 



and E-t representing the 'raising and lowering' operators of 5'C/_i(2). 
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3.2 Form factors: 0(2) versus SU-i{2) covariance 



Form factors in this context are matrix elements of some local quantum field between the 
vacuum and a multi-particle scattering state. They can in principle be computed from a 
recursive system of functional equations defined largely in terms of the given bootstrap 
S-matrix. Since the S- matrix has the enhanced 5'C/_i(2) symmetry it is natural to ask 
whether the associated functional equations are likewise covariant. Unlike the situation 
in other models this turns out to be not the case automatically, but it rather hinges on 
the values of certain statistics phases. Since this is a novel feature we briefly outhne the 
general structure of the form factor equations in the 0(2) bootstrap theory here. Details 
are relegated to Appendix A. Explicit results for some operators of interest are given in 
the next section and Appendices B and C. 

The form factors are tensors with respect to the obvious (real) action of 0(2) rotations. As 
with the S-matrix it convenient to diagonalize this action by the unitary transformation 
(3.5). We write /a„...ai(^n, ■ ■ ■ , ^i) for the components of some n-particle form factor 
in this "charged basis". The terminology is motivated by the fact that under a C/(l) 
transformation a form factor picks up a phase e^^'^, where e :— an + ■ ■ - + 0(1 plays the role 
of the U{1) charge. Equivalently e is the weight with respect to the Cartan subalgebra 
generator of SU -i{2). A form factor (of an operator) of Lorentz spin s should also have 
the homogeneity property 

/a„...ai (^„ + li, . . . , ^1 + k) = e^" (^n, . . . , ^l) . (3.8) 

For a fixed particle number n a form factor then has to satisfy the functional equations 

fan-aii^n, ■ ■ ■ , ^2, ^l) = Sl^^^^{(^2l) fan-ocsisi^n, ■ ■ ■ , ^1, O2) , (3.9a) 
/a„...ai(^n + StTZ, 6n-l, ■■■,02, 61) = F^^ fa„-i...ais(On-l, ■■■,&!, On) ■ (3.9b) 

In the second Eq. the shift by 27ri is to be understood in the sense of analytical continuation 
and is a constant matrix on whose role we elucidate below. First note that the system 
(3.9) decomposes into decoupled sectors with fixed U{1) charge e G {n,n — 2, ... , —n + 
2, —n} and dimension n!/n_!(n — n_)!, where n_ = (n — e)/2 is the number of '— ' labels in 
(«„, . . . , cti). Correspondingly the matrix F is diagonal in this basis but may be different 
in different charge sectors 

ri = Va{e)S^^. (3.10) 



11 



The phases i]a{e) can be thought of as statistics phases describing the relative statistics of 
the (quasilocal) operator whose form factors are considered and the field that generates 
the scattering states in a Haag-Ruelle type scattering theory [36]. See e.g. [29] for some 
simple examples. Iterating (3.9b) and employing the analyticity in u of (3.8) one finds 
the following spin-statistics relation 

r/_(e)"- r/+(e)'^-"- = e'"" . (3.11) 

A further condition arises if the underlying operator is hermitian. From fa„...ai {j^m ■ ■ ■ 
= f-ai...-a„{dl + ■ ■,&n + ^ne obtains 

Va{e)7j.a{-ey = 1 . (3.12) 



If only 0(2) invariance is assumed no further constraints exist and the phases rja{e) are 
part of the specification of a field operator in the bootstrap framework. Collectively they 
encode the information about the super selection structure of the theory. Since the first 
Eq. in (3.9) is covariant also with respect to the larger nonabelian SU-i{2) symmetry, it 
is natural to ask whether or not also the second Eq. is covariant for a suitable choice of 
the phases. The covariance requirement links the charge e sector with the e±2 sectors. It 
can be seen to entail an overdetermined set of relations for the phases i]a{e), - which turn 
out to be self-consistent. The requirement of quantum group covariance thus determines 
all phases ria{e) essentially uniquely; c.f. Appendix A. For n < 4 one finds explicitly: 

n = 2: r)+(2)^-ri±(0)^ri4-2) , 

n = 3: r/+(3) = TV±{1) = ±^±(-1) = ^-(-3) , 

n = 4 : 77+(4) = -r^±(2) = r^±(0) = -r^±(-2) = rj^-A) . (3.13) 

Generally, for fixed n, the relative signs are given by i]a{e) ~ exp y(e — na). Of course 
the actual phases solving (3.13) must be chosen n-independent. 

Let us illustrate the use of these relations in the charge e = 1 sector (where only the odd 
particle form factors are nonzero). We can take ^^±(1) = e'^^'^"' as the solution of (3.11). 
Then (3.12) fixes ti±{—1) = e^^'^^^. If we now in addition require that the field underlying 
these form factors is quantum group covariant, the e = ±1 sectors arc linked by (3.13), 
e.g. for n = 3. This yields the condition e^^'^*^ = —1, and we conclude: s = 1/4 mod 1/2. 
In words, an 0(2) doublet field that is in addition quantum group covariant can only have 
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Lorentz spin s = 1/4 mod 1/2. If we had started from the SU{2) invariant S-matrix (3.4) 
instead, no relative signs in (3.13) would have occurred, and an SU{2) doublet of 0(2) 
charge e = ±1 was forced to have Lorentz spin s = 1/2 mod 1/2, as expected. 

Next we proceed to the residue equations which link an n-particle form factor to an n— 2 
particle form factor. Consistency requires that the inverse of the matrix appears on the 
right hand side, irrespective of its concrete form. In the charged basis the precise formula is 
given in (A. 9). For generic phases (A. 9) will only be 0(2) covariant. Concretely this means 
that an n-particle form factor of [/(I) charge e is linked to an n—2 particle form factor with 
the same U{1) charge. For the choice (3.13) ensuring the SU^i{2) covariance of (3.9) one 
expects that also (A. 9) is quantum group covariant, which indeed turns out to be the case. 
Since all form factor Eqs then are covariant a quantum group transformation will map one 
solution into another solution, where "solution" actually means a sequence of functions 
whose members are linked by the (A. 9). Suitable sequences should correspond to (local) 
quantum fields in the 0(2) model. We thus find that field operators whose statistics 
phases enjoy the specific relation (3.13) form multiplets with respect to the quantum 
group action. Clearly one can concentrate on the multiplets transforming irreducibly; the 
multiplet and the associated form factor sequence will then be characterized by an isospin 
quantum number stemming from the representation theory of SU-i{2). In appendix A 
we list the irreducible multiplets for n < 4. 

We can summarize the situation as follows: The functional equations characterizing the 
form factors are 0(2) invariant and decompose into decoupled sectors with fixed U{1) 
charge e. In each sector statistics phases r]a{e) enter that are part of the specification of a 
field operator (or of an 0(2) multiplet thereof) in the bootstrap framework. In addition 
operators from different charge sectors whose statistics phases enjoy the particular relation 
(3.13) are members of an (irreducible) multiplet with respect to the nonabelian quantum 
group 5'C/_i(2). The existence of these multiplets is a nontrivial prediction of the bootstrap 
formulation. 

As remarked earlier there exists a (non-rigorous) transformation of the 0(2) model into 
a fcrmionic model with a manifest SU(2) invariance, for which the natural candidate 
S-matrix is (3.4), i.e. that of the SU(2) chiral Gross-Neveu model. As with the 0(2) 
- Sine-Gordon correspondence we expect that both systems share subsets of fields with 
identical correlation functions. An interesting one-to-one correspondence of the fields 
however is unlikely. To see this let us discuss the relation between the bootstrap systems 
based on the 5'C/_i(2) invariant S-matrix (3.1) and the SU{2) S-matrix (3.4) in more 
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detail: In the charged basis the mapping 

n 

l^n, ■ ■ ■ , ^l)a„...ai ^ |^n, ■ ■ ■ , ^l)a„...ai , (3-14) 

maps states whose exchange relations are governed by the SU (2) S- matrix (3.4) onto those 
whose exchange relations are governed by the SU -i{2) S-matrix (3.1), and also 'untwists' 
the non-trivial comultiplication of 5'C/_i(2) [43]. However (3.14) does not induce an inter- 
esting correspondence of the form factor sequences. For example if the SU (2) statistics 
phases are taken to be unity, the mapping (3.14) induces rf^^^^^{e) — e''^^^""")/^ for the 
SU-i{2) bootstrap system. This flips sign under n — > n — 2 while the statistics phases 
of a sequence acceptable for describing a SU-i{2) covariant field operator of course must 
be n-independent. Mathematically one can set up a correspondence between solutions of 
the form factor equations based on the SU{2) and the 5'C/_i(2) S-matrix. One way of do- 
ing this is to substitute the respective S-matrices into the Bethe Ansatz inspired integral 
formulae of [43, 3]. However this correspondence will in general not preserve the spin, 
the statistics phases, or even the covariance under the global symmetry group. One must 
conclude that there is no physically interesting one-to-one correspondence between the 
field content of the bootstrap systems based on the SU{2) and on the SU-i{2) S-matrix. 



3.3 Spin, Parafermion and Current form factors for n < 4 

With these preparations at hand we now seek to determine the form factors of the Noether 
current = -^{S^d^.S'^ - S'^d^S^) and the basic Spin field 5", a = 1,2. The latter is 
an 0(2) doublet and carries Lorentz spin s — 0. Prom the discussion following (3.13) we 
conclude that it cannot be a quantum group doublet as well. We thus also search for 
the form factors of an additional local "parafermion" field that is a, SU -i{2) doublet with 
Lorentz spin 1/4. (We shall comment on the relation to Smirnov's parafermion in the 
SG model [44] below.) Technically the construction of form factors for the Spin and the 
parafermion field is very similar. In order to treat both cases simultaneously we write 
^g{x) for the renormahzed field operators, with s — 0, 1/4, corresponding to the spin and 
the parafermion field, respectively. The objects of interest then are 

{O\J^{O)\9n,an;...;0i,ai) = -^e^ti. ^2!^p''(^i)j /«„... 01(6*71, 6*0 , neven,(3.15) 
(0|$^(0)|^„,a„;...;^i,ai) = C..«,(^n, • • • , ^1) , nodd. (3.16) 
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Here all components refer to the charged basis. The current form factors have charge 
e = while that of ^f{x) have charge e = ±1. The prefactor in the current form factor 
ensures current conservation; the on-shell momenta are po{0) = MchO, pi{6) = Msh^. As 
always in the form factor bootstrap Eqs (3.15), (3.16) must be regarded as a "statement 
of intent" . That is, the right hand is computed through the functional equation while the 
interpretation as the matrix elements aimed at on the left hand side has to be justified 
by additional considerations. 

As input for the recursive functional equations the normalization of the starting members 
has to be fixed. For the current a preferred normalization stems from the fact that the 
associated Noether charge Q should induce 0(2) rotations on the spins, i.e. [Q^-S"*] — 
itabS^. For the 2-particle form factor this converts into 

/+_(^2,^i) = -/-+(^2,^i) = ^ — + ..., e^^e^ + iT^, (3.17) 

6^2 — PI — ?7r 

where the dots denote regular terms. Writing f^^{62,9i) = /(^2 — ^i) the function f{6) 
has to satisfy the functional eqs f{9) = S2{9) f {~6) and f{9 + 2TTi) = —f{—9), with 5*2(6') 
from (3.1). They can be solved in terms of the function 

y{e) sh^e^W, 

r^cht(i +../.) -1 

^ ^ Jo t (1 + e*)sht ^ ^ 

which enjoys the following properties 

y{e) = S2{e)y{-e) , y{e + 2m) = y{-e) , 

y{9)y{9 + in) = - n_j_\y ' vi'^^) = ^ ■ (3-19) 

V2 27riJ \2TTi) 



We also note 



?7r 

A(0) = 0.304637 , lim [A(^) - A(-^)] = — . (3.20) 

0—*+oo 2 



Taking into account the residue condition (3.17) one obtains 



/+_(^2, ^1) = -f-402, e,) = ■ (3.21) 



With the 2-particle form factor explicitly known one can proceed to the 4-particle form 
factor. The formulas now get more involved and we defer the details to appendix B. The 
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Lorentz spin s = 1 is readily checked to be compatible with SU-i{2) covariance via (3.11) 
- (3.13). The (one-component) Noether current in the 0(2) model is the neutral member 
of a hidden quantum group triplet [43] , although the charge e = ±2 partners are nonlocal 
in the spin field. 

For the ^g{x) fields a natural normalization is 

(0|$^(0)|^,/3) = /^"(^) = 5fe-r (3.22) 

Proceeding to the 3-particle form factor we note that because of charge conservation and 
hermiticity there are only three independent components 

fl{Ol, 92,03) = fl-+{03,02,0i) , 

f3{Ol,02,es)^f-__{03,02,ei). (3.23) 

The general 3-particle residue equation (A. 9) in the charge e — ±1 sectors 

^res,3,,=,./,"3,^,^(^) = ^«3+7 hie)'' 542«,(^2i) - fpi0i) , (3-24) 

thus translates into 

fk{u,v,0)^ '^^—Wk{v-0), u^v + iir, (3.25) 

U — V — ZTT 



/++_(^3, 02, 0l) — 
U-+{03,02,0l) = 
f-++{03, 02, 0l) = 



where with r] := r7+(l) = e^'^'* and the notation from (3.6) 



Wi{9) = -iiSR{e) , W2{0) = I - iiSt{0) , W3{0) = i-r^S2{0) . (3.26) 

These functional equations can be solved for generic s, the solution is described in ap- 
pendix C. For the spin field s — Q fixes our candidate form factors. For the parafermion 
field we know that SU -i{2) covariance (regardless of irreducibility) requires s = 1/4. (The 
specific construction used in appendix C removes the additive mod 1/2 ambiguities). As 
explained in appendix A the condition that the parafermion field (and hence its form 
factors) transform irreducibly as a isospin 1/2 doublet requires in addition 

C(^3, 02, 0l) := fl{03, 02, 0l) - f2{03, 02, 0i) + ^(^S, 02, ^l) = . (3.27) 
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For the solution constructed in Appendix C, C(^3) ^2, ^1) can be shown to be proportional 
to rf + 1, so that s = 1/4 also entails the desired irreducible transformation law. 

So far we didn't say anything about the local structure of the parafermion field ^"^^(a;) 
supposed to underly the above solution of the form factor equations. We can address 
this point by employing a result by Smirnov [44] on the existence of parafermionic fields 
in the Sine-Gordon model. Smirnov's fields have well defined exchange relations of the 
form ^a{.x)^h{.y) = 'S'ab(±C)o)^(;(t/)\E'c(x), for ±.x^ > ±|/\ where x^^y^ are the space 
components of y in a fixed Lorentz frame. For generic Sine-Gordon coupling (3 these 
fields are nonlocal because the two limiting S-matrices are distinct. From the analysis of 
the form factors of the energy momentum tensor one can see that our parafermion field 
is the Stt limit of Smirnov's field. But thanks to (3.2) it is now a /oca/ field. Indeed 

for the charged components the exchange relations assume the simple form 

*^/4(^)^V4(^) = -%A{yM/S^) , (3-28) 

for all spacelike separated points y. Likewise the transformation properties under the 
quantum group change qualitatively. As analyzed by several authors [44, 30, 16] in general 
a dynamical quantum group symmetry of the S-matrix acts in a nonlocal way on the 
field operators. In contrast the field ^"^^(x) transforms nicely as an irreducible 5'C/_i(2) 
doublet. The situation is thus reminiscent of the Ising model, which can be viewed as 
the n = 1 case of (2.1). There both the spin and the fermion are local fields of Lorentz 
spin and 1/2, respectively. Both are relatively non-local but generate equivalent sets 
of scattering states. Though in the absence of a field theoretical construction of the 
parafermion field it is difficult to examine this point, we expect the interplay between the 
spin and the parafermion field in the XY QFT to be analogous. 



17 



4. Lattice computations 



4.1 General setup 



For the lattice regularization we consider a square lattice with action 



S=-Kj2S{x)-S{x + fi), 



(4.1) 



where S{x) ■ S{x) — X^a '^"(^)'^"(^) ~ ^- Solving the constraint with S^{x) — cos(p{x), 
S'^{x) — smip{x) the action reduces to the standard XY lattice action Eq. (2.18). The cor- 
relation functions are defined as in the continuum theory except that the spatial integrals 
are replaced by discrete sums. There is an enormous literature, both on numerical sim- 
ulations of the XY model [19] and on its high temperature expansion [10, 40]. Presently 
the best numerical estimate of the critical point for the standard action is [21] 



Previous numerical investigations mostly concentrated on the comparison with the KT 
theory. We will outline the aspects relevant here and some refinements in Section 6.1. 
Our main goal however is to compare the continuum limit of the XY model with the 0(2) 
bootstrap theory. 

In the rest of this section we collect some general information on our simulations and 
continue with detailed discussions of the measurement of various observables in Sections 
5 and 6. All numerical simulations were done on an L x T lattice with periodic boundary 
conditions in each direction. During the entire investigation two independently written 
programs were employed and many cross checks were made. Both used multi-cluster 
updating. The Ising spins are embedded hke in Wolff's single cluster algorithm [50]; the 
resulting Ising model is then updated with a generahzation of the Swendson-Wang [46] 
multi-cluster algorithm. In one application of the program each run started from a random 
configuration and consisted of a large number of sweeps of which a large initial proportion 
were used solely for equilibration. In another application after initial equilibration the final 
configuration of the run was stored and read in for the starting configuration of the next 
run. In most cases the observables were measured using improved (cluster) estimators. 
The final data sample of the many runs was averaged and the error was computed using 
the jack-knife method. 



= 1.1199(1) . 



(4.2) 
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Since we aim at achieving high precision, for many quantities to an accuracy of < 1%, a 
considerable source of concern to us was the random number generator (RNG). Indeed 
at an initial stage of this project we found that results obtained by various RNGs could 
differ by many standard deviations. We thus subjected the RNGs to several tests, specific 
to the model and the quantities considered here. The results are reported in Appendix E. 

Numerical simulations of course are restricted to work in finite volume. In the next two 
subsections we therefore discuss finite volume effects, first in the continuum and then on 
the lattice. 



4.2 Finite volume effects in the continuum 

If a continuum QFT is confined in a box the physical observables will depend on the 
geometry and boundary conditions. Consider first the mass m(L) of the 1-particle state 
on a cylinder of circumference L. Liischer [32] has shown that in a theory without a 
"three-point coupling" of this particle to itself or to any other single-particle state, for 
large physical volumes z — ML ^ 1 the finite volume dependence of the mass is given 
by: 

= !!<L)-A£ ^ 1 r ^-.cosht cosh(i)/(i) + . . . , (4.3) 

M vr Jo 

where f(t) is the forward scattering amplitude at an off-shell point. Thus in these models 
the infinite volume limit is approached rapidly at a rate ~ exp{—z) (from above if 
/(O) > 0). 

More explicitly in the 0(n) sigma-models the amplitude / is given by 

f(t)^n-^\2So{9)+n{n-l)Si{9) + {n + 2){n- 1)32(9)] . (4.4) 



2n 



=t+iTc/2 



With the proposed S-matrix (2.12), (3.1) of the 0(2) model one obtains 



9-2 

m = ~ - 



^2 + 7^2/4 

in which case Eq. (4.3) results in 



r(3/l-/Y/2^) 



r(l/4-ii/27r) 



(4.5) 



D{z) ~ 1.162475182 ^ [l + 0.36432/z + 0{l/z'^)] . (4.6) 
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Next we consider the zero-momentum coupling in a square box of length L in each direc- 
tion. To get a feeling of the finite volume effects consider the expression for gji{L) in the 
leading order 1/n expansion [12]. One obtains 

g^{L) = ^r(oo)[1 - V8^z^^e~%l + 0(l/z)) + •••], (4.7) 

with A4 = 1/2. Although the 1/n expansion is not expected to be quantitatively applica- 
ble to n = 2, one might hope that the qualitative features are correct, in particular that 
the finite volume effects are exponentially suppressed and secondly that the exponent /I4 
of the multiplicative power correction is independent of n. Some corroboration of this 
might come from investigating the finite volume effects in a square box in an effective 
Lagrangean framework similar to [32]. 



4.3 Finite volume effects in the lattice regularization 

The particular lattice sizes used will be specified later, but they were generally selected 
to enable studies of finite size effects at fixed correlation length and vice versa, subject of 
course to restrictions due to the CPU power available to us. 

For a coupling gR{K, L) depending on the bare coupling K and size L x L which tends 
to gji in the continuum and infinite volume limits, we adopted the 0(n) definition 

,K(i^,L)=f^Vh + 2-lgX^l, (4.8) 



for n = 2, where = S"'{x). In Eq. (4.8) and throughout Sect. 6, ^ is taken to be 



' Z'^C'-l, (4.9) 



2sin(7r/L) y ^(fco) 

where ko = (27r/L, 0); c.f. ref. [15]. This correlation length converges in the thermody- 
namic limit to the second moment correlation length 1/Mr. In this connection we would 
like to draw the reader's attention to a subtlety which is discussed at the end of Section 
4 in ref. [7] . 

The Noether current on the lattice is defined by 

J^{x) = K {S\x)d^S\x) - S\x)d^S\x)] , (4.10) 
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where d^f^x) = f{x + fi) — f{x). Introducing its two point function as 



(4.11) 



X 



with q — {qi, ^2), the well-known Ward identity (for the standard action) reads 



(4.12) 



where the energy expectation E is 



E{K,L) = Y,{Six)-S{x + fi)). 



(4.13) 



Next we wish to define for a finite periodic lattice the counterpart of the transverse 
current correlation function /(A;^) in (2.11). J^j, can naturally be decomposed into 3 
pieces ('transverse', 'longitudinal' and 'harmonic'), as discussed in [39]. The harmonic 
piece is concentrated at the origin in momentum space and has the value Jii(0|A', L). 

In the thermodynamic limit, because of the presence of a mass gap, J^y{q\K^L) will be 
a real analytic function; so it cannot contain any remnant of the harmonic piece (which 
would be proportional to a 5 function). The longitudinal and the transverse parts have 
to go to the same limit as g ^ to avoid any non-analyticity at zero momentum. Being 
a contact term the value of J^;y(0|-ft', L) has no physical meaning. For this reason, as 
explained already in Section 2.1, we define the transverse part in such a way that it 
vanishes at zero momentum in the thermodynamic limit. This suggests the definition 



to which we shall refer as the SUB definition. It ensures /(O) = at finite L and ^. 

Another possible definition is suggested by the Ward identity, which for momenta of the 
form q = (0, 52) reads 



We can use (4.15) as an alternative definition for a lattice counterpart of /(A;^) in (2.11), 
to which we shall refer to as the WARD definition. The normalization /(O) = is then 



/((O, q2)\K, L) := Jn((0, Q)\K, L) - Jn((0, q2)\K, L) , 



(4.14) 



/((O, q2)\K, L) -E{K, L) - Jn((0, ^2) L) . 



(4.15) 
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only restored in the thermodynamic hmit, but for numerical purposes the WARD version 
often is advantageous. To verify that (4.15) vanishes at g = in the thermodynamic limit 
note that Eq. (4.12) also entails 

Ju{{qi,0)\K,L) = ^EiK,L), Vgi ^0. (4.16) 

Since Jii{{qi,0)\K, L) becomes a real analytic function of gi e [— tt, tt] for L — > oo 
Eq. (4.16) remains vahd also for qi — [39]. In the thermodynamic hmit therefore 
the WARD and the SUB definitions are equivalent. In Section 6 we shall consider both 
options. 

Unfortunately few rigorous results exist on finite size effects ior z — L/^{K, L) ^ 1 on 
the lattice. In general we expect in the 0(n) models, for a large class of observables O 
and for fixed bare couphng K, the finite size effects to be either of the form 



-z 



Co^(X) + Cf(X)- + 0(lA^) 
z 



+ ..., (4.17) 



or the same with Cf{K) replaced by 0[K, oo)Cf{K). Here either the amplitudes Cf{K) 
or OiK^ oo)Cf{K) are hoped to be almost constant for large correlation lengths. Further 
the dynamical assumption enters that the fall-off is e~"^ with a = 1, as expected in the 
absence of two-particle bound states. Further it seems reasonable to expect A'-^ to depend 
only on the form of the observable and not on the dynamics e.g. not on n. 

As mentioned earlier one can probably, for certain correlation functions, better justify 
these assumptions by performing a Feynman diagram analysis in the framework of an 
effective massive lattice field theory analogous Liischcr's in the continuum theory [32]. 
Here we only indicate some plausibility considerations based on the leading order in the 
lattice 1/n expansion (viewed as a summation of bubble diagrams) . For the coupling this 
gives 

ng^{K,L)^2A{0)/ml + Oil/n), 

with A(0)-^ = -i5^(E, + m^)-^ i^ = ;^E(^P + ^o)"'- (4.18) 

p p 
Here ^ = l/rrio and the sums range over momenta = 27in^/L, with < < L — 1, 
11 — 1,2, and Ep = 4 sin^p^/2. For the current two point function the leading term is 





e — e' 




^Ep + ml 


Ep + ml 
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where P = p + q, and for the energy expectation it is 



E{K,L) 



k1? 



E 



p 



E^cosp^ 



(4.20) 



Starting from these formulae one can study separately the finite volume effects and lattice 
artifacts in this approximation; the results confirm the before mentioned assumptions. We 
stress again the big qualitative difference between the finite volume effects and the lattice 
artifacts. Whereas the finite volume effects represent continuum physics and hence are 
expected to be structurally universal, the lattice artifacts are in general non- universal. In 
particular for the spin 2-point function one finds in this framework = = — 1/2 for 
the exponent in (4.17), whereas for the current and g-^i we get A'^ = ^4 = 1/2. (Both the 
current 2-point function and g^i depend linearly on spin 4-point functions and perhaps 
this accounts for the same exponent A^.) If one does not wish to adopt this framework 
the exponents can be kept as fit parameters; c.f. Section 6. 

4.4 Lattice eirtifacts 

After the extrapolation to infinite volume has been performed, the results can be regarded 
as corresponding to a lattice 0(2) action in infinite volume. The extrapolation to infinite 
correlation length is usually hampered by the lack of information about the rate of ap- 
proach. Based on the Sine-Gordon description of the KT transition [1] one of us [4] has 
argued that for the XY model, say with standard action (4.1), the leading lattice artifacts 
are calculable from the continuum Sine-Gordon theory. This applies to obscrvablcs like 
the S-matrix or the two-point function of the Nocthcr current, where already at finite 
correlation length a preferred normalization exists. Implicit in this proposal is a certain 
degree of action-independence of the leading lattice artifacts, but at present it is not clear 
to which class of actions it applies. Later on we test this proposal for the standard action 
and the two observables mentioned. 

Let Uxy{C) denote such an observable, where ^ is the correlation length in lattice units 
and the dependence on other variables is suppressed. Then t/xY(C), computed e.g. with 
the standard action is predicted to be of the form [4] 



Here uq is the continuum value and Ui is the leading correction. The parameter u is 
action-dependent but should not depend on the physical quantity considered. Note the 




+ O((ln0-^). 



(4.21) 
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extremely slow decay ~ l/(ln,^)^. As explained in [4] the peculiar structure of the KT 
phase diagram allows one to relate both Uq and Ui to the continuum SG theory. Namely 
if Usg{'^) denotes the counterpart of the physical quantity considered in the continuum 
SG theory with coupling u close to u = 0, then 

Usg{^) ^uo + uiu^ + 0{u^) . (4.22) 

The equality of the leading Uq term in (4.21) and (4.22) simply reexpresses the link 
between the XY model and the SG QFT alluded to in the introduction. Remarkably the 
coefficient Mi of the first correction is likewise the same in both cases. 

Our first application of these formulae is to the scattering phase shifts. Recall the SG 
model S-matrix (2.12) with (2.16). Introducing the phase shifts by ^'/(^li/) = ex.p{2iSi{9\i')), 
1 — 0,1,2, and expanding in i/ at fixed 9 yields 

6i{e\iy) = 5i{e) + vH'iie) + o{p^) . (4.23) 

By construction 5i{9) are the phase shifts of the proposed 0(2) S-matrix (3.1). The 
0{v'^) coefficients can be related to the lattice artifacts by the relations (4.22), (4.21). 
They come out as 

This is used to predict the leading lattice artifacts in the phase shift analysis of Section 
5. 

It is also feasible to calculate the leading lattice artifacts for the two-point function of the 
Noether current. In [4] this was done in two-loop perturbation theory. Alternatively, using 
the current form factors of the SG model, the leading artifacts can also be calculated non- 
perturbatively via the form factor bootstrap. We worked out the two-particle contribution 
to this correction. The comparison with numerical data is presented in Section 6. 
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5. MC results for the phase shifts 



We begin by numerically investigating the S-matrix. It is the prime input for the boot- 
strap formulation and an appreciable discrepancy to the bootstrap result (3.1) would 
immediately rule out that the 0(2) bootstrap theory describes the continuum limit of 
the XY model. The technique to numerically determine the S-matrix takes advantage 
of the fact that the large volume dependence of the spectrum in a periodic (spatial box) 
encodes information on the infinite volume S-matrix. For example the volume dependence 
of the (stable) 1-particle mass is governed by the forward scattering amplitude [32], and 
a determination of the low-lying two-particle spectrum gives a measurement of the low 
energy two-particle phase shifts [31]. Here we restrict the discussion to a field theory in 
1-1-1 dimensions, i.e. the "spatial box" in this case is just a circle of circumference L. To 
our knowledge the first attempt to determine the phase shifts of the XY model was by 
Vohwinkel [49], and we record his results in Appendix F. 

5.1 1-peirticle mcisses 

We chose to measure on the same lattices as Vohwinkel, firstly since they are practical and 
secondly for having the advantage of being able to compare independent measurements. 
These lattices arc listed in Table 1 together with the measured 1-particle masses. In each 
case the "time" extent of the lattice is T = 2L and periodic boundary conditions are 
imposed in each direction. 



K 


L 


m(L) [49] 


m{L) 


0.86 


64 


0.1711(1) 


0.17096(4) 


0.92 


128 


0.09465(3) 


0.09461(6) 


0.97 


256 


0.04620(1) 


0.04603(14) 



Table 1: Values of K and L used in the measurements, with the 1-particle mass in{L) [49] obtained by 
Vohwinkel [49] and our measurements, m{L). 

What is quoted are the results for the single particle masses obtained by fitting the zero- 
momentum spin 2-point function with a 2-mass formula, with the second mass constrained 
to be 7712 = 3mi. (The 1-mass fit and the unconstrained 2-mass fit give the same values for 
mi within the errors.) The three different values of K correspond to correlation lengths 
~ 6, 11, 22. In all cases there is good agreement with the results of Vohwinkel. 
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All the lattices are chosen to have m{L)L > 10, so that finite volume effects on the 
1-particle masses are expected to be very small according to Eq. (4.6). We denote the 
resulting 1-particle mass by m. 



K 


L 


m(L) 


D 


-^theor 


0.86 


32 


0.17135(3) 


0.0023(4) 


0.0022 


0.92 


64 


0.09478(3) 


0.0018(10) 


0.0012 


0.97 


128 


0.04622(4) 


0.0041(39) 


0.0014 



Table 2: 1-particle masses and finite volume effects. 



In order to quantitatively test the latter expectation we measured the 1-particle masses on 
lattices with the same three bare couphngs K as those in Table 1 but with half the previous 

spatial extent. Our results are listed in Table 2. Also shown are the measured values of 
the finite volume mass shifts D = {m{L) — m{2L))/m{2L) and -Dtheor computed from 
Eq. (4.6). At the smallest correlation length the measured shift is completely consistent 
with our ansatz for the candidate S-matrix. However, lattice artifacts are to be expected 
and unfortunately our measurements at the larger correlation lengths are too imprecise 
to study these effects. 

5.2 Phase shifts 

For a numerical test of the proposed S-matrix (3.1) it is useful not to presuppose the 
symmetry enhancement. That is we adopt the generic 0(2) invariant parameterization 
(2.12) with n — 2. In terms of the Si{9) the phase shifts are defined as 

Si{e)^ew{2i5i{e)] , (5.1) 

and can simplified to 

S2(0) 

The last relation is responsible for 



TT U 

= 6i{9)-\ arctan — , (5.2a) 

2 TT 

= 5,(0) -|. (5.2c) 

the symmetry enhancement discussed in section 3.1. 
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To measure the phase shifts on the lattice one sets out to determine the center-of-mass mo- 
menta of 2-particle eigenstates of the transfer matrix, since due to the periodic boundary 
condition in the spatial direction these momenta are quantized according to 

PnL + 25 (On) ^27771, p„ = msinh^^„. (5.3) 
To accomplish this one measures correlators 

Cg\t) = {Oj{x,0)Oj{y,t)), (5.4) 

where the O^^^ are 2-spin operators with zero total momentum in the isospin channels 
7 = 0,1,2: 

1 

Oi{x, t)^- t)S''i^ + x,t), 6o = 1 , &1 = &2 = 2 , (5.5) 

z=0 

the Pi being the projectors in Eq. (2.14). Taking T large enough so that terms propor- 
tional to e~^"*^ can be neglected, one has 

n 

where 

iji!\x) = {vac\Oi{x,0)\n) , (5.7) 

is the "wave function" of the corresponding state. In the / = channel the vacuum 
also contributes as an intermediate state. In this case one can subtract this contribution 
from the beginning, i.e. consider the connected correlator. Alternatively, one can take the 
full 7 = correlator, keeping in mind that in this case the vacuum is the lowest energy 
intermediate state. 

Now there are at least two ways to proceed. The first is to extract the energies En of the 
2-particle states which dominate the correlation function Eq. (5.6) for sufficiently large t, 
and then compute the corresponding center of mass momenta via 

Er, = 2^(1) (p„) = 2Vp2 +m2 . (5.8) 

This was the strategy used in the pioneering paper of Liischer and Wolff [31] and adopted 
by Vohwinkel in his studies [49]. In Eq. (5.8) lattice artifacts have been neglected and the 
physical volume is taken so large that the finite volume dependence of the single particle 
masses is negligible. 
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In ref. [6] an alternative was suggested which starts from the observation that the relative 
momentum 2p„ of the two particles is also encoded in the wave function: in the symmetric 
channels (/ = 0,2) one should have 

il^nix) = Acospn{x - L/2), foT R < X < L - R , (5.9) 

and similarly with sinp„(x — L/2) for the 1=1 channel. Here R is the "interaction 
range" characterized by the requirement that for a relative distance x > R the two 
particles propagate essentially freely. Note that Eq. (5.3) assumes that the box is large 
enough to accommodate the two particles without "squeezing" them, i.e. L/2 > R. 

The rank N of the matrix C{t) in Eq. (5.6) is L/2, L/2 - 1 and L/2 + 1 in the / = 0, 1, 2 
channels, respectively. (This is when the connected correlation function is considered in 
the / = channel, otherwise N = L/2 + 1 also in this channel.) We assume that for t > to 
(with some to) no more than N states contribute to CxJ{t), i.e. that the contribution 
from the states n > N can be neglected completely. 

Liischer and Wolff [31] suggested to determine the energies En from the generalized eigen- 
value problem* 

C{t)Vn^Xn{t,to)C{to)Vn. (5.10) 

Provided the sum in Eq. (5.6) is restricted to N terms, 1 < n < N, the eigenvalues of 
Eq. (5.10) are given exactly by 

A„(t,to)=e-^"(*-*°). (5.11) 
It is easy to show that (apart from the normalization) 

V'n(a^) =^aj/(to)^^n(y)- (5.12) 

y 

A problematic feature of the generalized eigenvalue equation (5.10) is that its solutions 
become unstable if C{to) has very small eigenvalues. This can be seen by observing 
that A„(t,to) are the eigenvalues of the ordinary eigenvalue equation for the matrix 
C{tQ)^^^'^C{t)C{tQ)^^^'^. Of course, the exact correlation matrix C(to) is positive defi- 
nite, but the statistical noise will spoil this property, and the measured matrix can have 
even negative eigenvalues, especially for larger values of to and for large number of opera- 
tors A^. For this reason in ref. [31] ~ L/4 operators were used (actually, in momentum 
space rather then in x-space) and the values of to were restricted to and 1. 

*This equation was considered already earlier by Michael [35], in connection with a variational ap- 
proach evaluating the static potential in lattice gauge theory. 
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To avoid the instability, we restrict first the correlation matrix to an M-dimensional sub- 
space (M < A^) spanned by the first M eigenvectors of C(to) with the largest eigenvalues 
(still stable against the statistical fluctuations) [6, 37]. The generalized eigenvalue prob- 
lem is then written for the new correlation matrix C{t) in this reduced basis. Of course, 
to read off the momenta, the wave functions have to be transformed back into the original 
basis labeled by the relative distance x. 

Following refs. [35, 31] one can obtain En from the plateau of the "effective energy" 

(t) = In ^ ^ , (5.13) 

A„(i + l,to) ^ ' 

and determine the corresponding momentum p„ from Eq. (5.8). 

The alternative way is to fit the wave function ipn{x) by the ansatz in Eq. (5.9) for 
< X < L/2. We have verified that it is safe to take Xq > 3/m. There is also a large 
window where the results are not sensitive to the variation of M, for different choices of 
to. For the largest correlation length ,^exp ~ 21.6 we could take to as large as 10, which 
would be impossible without the preceding truncation. 

In [6] we concluded that the wave function method has somewhat smaller errors and is 
more stable. In particular the smallest momentum obtained from the 2-particle energy is 
quite sensitive to the error in the single particle mass, while in the wave function method 
the value of this mass is not used at all. Although we measured the phase shifts by both 
methods and checked their consistency, only our results from the wave function method 
will be presented here. Vohwinkel's results on the energy levels are recorded in Table 15. 

Fig. 1 shows the wave functions of the first six 2-particle states in the / = channel 
obtained using Eqs. (5.10), (5.12). Fig. 2 displays the deviations of the first three wave 
functions from the corresponding free one, Acosp{x — L/2) in the 1 = channel. As 
expected, the true wave functions deviate from the free ones only for small relative dis- 
tances of O(^cxp)- The plots illustrate that the momentum p can be determined quite 
precisely by fitting the wave function in some properly chosen range xq < x < L/2. Note, 
however, that in Eq. (5.3) the momentum pn is multiplied by L, hence it has to be de- 
termined to good precision in order to yield a reasonable error for the phase shift. In our 
simulations therefore only the phase shifts at the first 3-4 momenta could be determined 
with a reasonable error. 
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Figure 1: The first six wave functions (vac|0/(x)|n) in the J = channel for K = 0.97, L = 256. 
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Figure 2: Deviations of the first 3 wave functions from the ansatz Acosp{x — L/2). 



A summary of our measurements of the phase shifts is given in Tabic 3. The Figures 
3 to 5 compare these results with the theoretical curves of Eqs. (5.2a-c). The leading 
lattice artifacts according to Eq. (4.24) for the largest correlation length ^cxp = 21.645(5) 
are also shown. The present overall results are certainly consistent with the theoretical 
expectations. 
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Table 3: The phase shifts 6i{p) for the first 3 states (n = 1, 2, 3) in different isospin channels, with the 
values of p/m determined from the wave function. The continuum results [p/m]ex and [^/(p)]ex calculated 
from the S-matrix are also given. 
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Figure 3: 7 = phase shifts. Points denoted by circle, square and diamond correspond to correlation 
length w 6, 11 and 22, respectively. The solid line is the continuum result, Eq. (5.2a). 



0.4 




Figure 4: / = 1 phase shifts. The notations for data points are the same as in Fig. 3. The solid line is 
the continuum result, the other two lines include leading corrections due to finite ^exp of Eq. (4.24) with 
u = 1.46 (see [21]): the dashed line for ^exp ~ 22, the dot-dashed line for ^exp ~ H- 
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Figure 5:7 = 2 phase shifts. The notations are the same as in Fig. 4. 



33 



6. MC results for correlation functions and 4-point coupling 

Here we describe our MC results for the intrinsic 4-point coupling gji and for the current 
and spin 2-point correlation functions in Fourier space. In addition we reconsider some 
aspects of the Kosterlitz-Thouless (KT) theory. 

A list of the lattices considered together with some of the MC results is given in Table 4. 
We group the lattices into famihes labeled 1 to 12. Members of a given family correspond 
to the same coupling K but different size L. Throughout this section ^ denotes the 
second moment correlation length (4.9), which is a function of K and L. We denote the 
'apparent' physical size of the lattice by z' = L/^{K, L). For most lattices we made 200k 
measurements where we measured all the physical quantities above supplemented by an 
additional 2M measurements for the "bulk" quantities (,^, x ^"^^ 9b) only. We also took 
data on 'thermodynamic' lattices {L/^ ^ 14). These are reported in Appendix E. We 
did not use them in our analysis because there were some not completely understood 
inconsistencies between data taken on different machines, with different random number 
generators and different programs. 



6.1 KT theory 

We begin by studying the dependence of the correlation length ^ and the susceptibility 
X on the coupling K and compare our MC results to the predictions of the KT theory. 
For this purpose we first need to extrapolate the measured values of ^ and x to infinite 
volume. This is done by a finite size scaling analysis. 

For fixed coupling the size dependence of the correlation length is assumed to be given by 



Cl{K) + Ci{K)l 
z 



+ ... 



(6.1) 



Here z := L/^cxp{K,oo) where ,^cxp is the true ("exponential") correlation length, ^exp 
is expected to be very close to the second moment correlation length: the form factor 
bootstrap predicts ^cxp/^ = \fyif^2 = 1.00089 and the lattice models are certainly not 
very far from the form factor construction; therefore at our ^ 0.001 accuracy the distinc- 
tion between the two correlation lengths can be neglected. For the exponent in (6.1) we 
expect = ^2 = — 1/2 on account of the considerations outlined in Section 4.3, but wc 
also tried different values of that exponent in order to see if the data indeed support this 
expectation. 
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Table 4: Lattice parameters and results. 
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We produce a global fit of the finite size effects by fitting all the data in Table 4 to the 
form (6.1) truncating after the Cg term and taking Cq to be independent of K. The 
values of ^{K, oo) are fit parameters, one for each of the twelve families. Since initially 
we do not know the values of L/^{K, oo) we first replace it by L/^{K,L); this leads to 
a first estimate of ^{K, oo) which is then used in the fit ansatz. Iterating this procedure 
about 5 times leads to convergence of our extrapolated values of ^{K, oo). It turns out 
that this type of fit favors a value of A2 near —1/2 in agreement with expectation; for 
this value of the exponent we obtain chi^/dof = 0.9 with 19 degrees of freedom. 

Furthermore, including a term makes the fit very insensitive to the value of the ex- 
ponent A2, both as far as the fit quality and the extrapolated values of ^{K,oo) are 
concerned. For A2 = —1/2 the coefficient comes out consistent with zero. Therefore 
in Table 5 we only report the values obtained with the simplest fit with A2 = —1/2 and 
only one finite size correction term. 

Figure 6 illustrates the z dependence of the correlation length for K — 0.97. 

21.7 
21.6 
21.5 
21.4 
21.3 
21.2 
21.1 

0.002 0.004 0.006 0.008 0.01 

Figure 6: ^{K,L) vs. exp{-z)/^/z at K = 0.97. 

The FS analysis of the susceptibility is analogous, but we don't need any iteration, since 
we start already with the right z values. Again the data favor a value near —1/2 if we 
truncate with the Cq term, and if we include the next term, the fit becomes insensitive 
to A2. In Table 5 we present the results from the simplest fit assuming A2 = —1/2, which 
has a chi^/dof of 0.7 with 19 degrees of freedom. 

Next we study the A'-dcpcndence of the infinite volume quantities. One of the best 
known predictions of KT theory [28] is the unusual coupling dependence of the correlation 
length. Close to the critical point Kc, ^ is predicted to diverge as ^ ~ exp(6/i/r), where 
T = Kc — K is the reduced coupling and 6 is a non-universal constant. In more detail, 
the Sine-Gordon description of the KT transition entails (using eqs. (8.13) and (8.14) of 
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[1] *) 

\n^{K,oo) = -^-u + c^/¥+ (6.2) 

where u and c are again non-universal constants and the dots stand for higher powers 
in T. We fitted the ^{K,oo) data to the form (6.2) without higher terms, leaving out 
different numbers rigUp of the low ^ families to see how stable the fit parameters are. The 
results are presented in Table 6. Note that some of the fits have an unacceptable chi^. In 
any case, the determination of u is not very stable. This situation could be a sign of the 
inappropriateness of the ansatz (6.2), or of some problem with our data or it could mean 
that we are still too far from the critical point, so that asymptotic formulae cannot yet 
be applied reliably. In further fits below that do involve u, we use the value appropriate 
to the number of discarded families. A visual illustration of the ngkip = 1 fit is shown in 
the left part of Fig. 7. 

We should mention that Hasenbusch and Finn [21] have determined the constants b and 
u by their method of matching to the exactly solvable BCSOS model; their values are 
u — 1.46(1) and b — 1.879(4), in rough agreement at least with some of our estimates. It 
should be noted, however, that their method avoids the problem of controlling subleading 
corrections like the third term in eq. (6.2). 
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Table 5: Correlation length and susceptibility extrapolated to infinite volume using a global fit. 



KT theory also predicts the asymptotics of xiO large ^ (at least up to possible loga- 
rithmic corrections), 

X^e~\ with 7^ = 1/4. (6.3) 

A hnear fit of Inx versus In^ is shown in Fig. 7. The slope is 1.73, which is very close to 
the expected result. The fit is visually good, but even if we omit famihes 1,2 and 3 we 
get a huge value of chi^/dof ~ 288/7, indicating the presence of non-negligible subleading 
terms. Irving and Kenna [25], following Butera and Comi [11] (see also [22]) argued that 

*note that (8.15) contains a misprint. 
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Table 6: Fit of the infinite volume correlation length to the ansatz (6.2); rigkip denotes the number of 
low ^ families discarded. 
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Figure 7: Illustration of KT scenario: Left, In^ versus r and a fit to (6.2). Right, Inx versus In^ and 
a linear fit. 

the Kosterlitz-Thouless theory implies the following refinement of (6.3): 

X-e-\\^iV\ (6.4) 

with r = —1/16. This would in particular mean that x grows faster than ^^'^^j whereas 
the data on the contrary indicate a slower increase. 

On the other hand one of us [4] has argued that the Kosterlitz-Thouless theory imphes 
r = 0, with a specific additive correction to (6.3): the renormalization group invariant 
quantity Q was introduced, which for large ^ behaves as 
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It was then argued that the correct asymptotic formula, instead of (6.4), is 

lnx~^lne + 0(Q). (6.6) 

Taking u from Tabic 6 the relation (6.6) can again be tested against the data. We fitted 
In;^ — 1.75 In,^ against Q, discarding successively more and more low ^ families; in the fits 
we used the best value of u corresponding to the same number of discarded families. The 
results are given in Table 7. 
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44/9 


32/8 


8.2/7 


5.4/6 


2.9/5 


0.6/4 


0.5/3 


0.5/2 



Table 7: Fits of In x — 1.75 In^ to linear function of Q; rigkip denotes the number of low ^ families 
discarded. 

Starting with rigkip = 3 the fits are acceptable, but of course it should be remembered that 
we are not quite sure which value of u should be used. We do not list the fit parameters, 
but they are quite stable. So one can say the data support the prediction (6.6). 
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Figure 8: Inx — 1.75 In ^ versus Q and a linear fit on the data from the ^ > 39 lattices. 

We will use the variable Q as the quantity characterizing the distance from the continuum 
limit. 

6.2 Determination of 

Lattice determinations of (?r can be based cither on the high temperature expansion or 
on numerical simulations. The results obtained via the high temperature expansion and 
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the standard action (4.1) are 

= 9.15(10) [10] = 9.10(5) [40] 



(6.7) 



In the numerical simulations we again aimed at achieving a precision of better than one 
percent. For the necessary extrapolation to the infinite volume we use the procedure 
outlined before, i.e. we fit the data to the ansatz 



In gn{K, L) = \ng^{K, oo) + z^^e 



Co^'^(X) + Cf (X)^ + ... 



(6.8) 



where, as described in Section 4.3, there are arguments suggesting A^^ = = 1/2. As 
opposed to the situation with ^ and with this value, the leading finite size correction 
alone docs not properly describe the FS dependence. If we instead use the optimized 
value A4 = 0.8, a subleading finite size correction is not needed. It is gratifying that the 
extrapolated values are almost independent of which of the two options we choose. We 
report the infinite volume values obtained with = 1/2 and two FS correction terms as 
well as those with only the leading FS correction and A4 = 0.8 in Table 8. Both fits have 
a chi^/dof around 1. 
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Table 8: Fits to with A4 = 1/2 using 2 finite size parameters (upper line) and with = 0.8 using 
only the leading finite size parameter (lower line). 

Finally we turn to the continuum limit of g^i- One of us [4] has argued that the leading 
lattice artifacts are proportional to the quantity Q, which in turn depends on the param- 
eter u extracted from the fits in table 6. We present in Table 9 the results of fits of the 
infinite volume g^ values to a linear function of Q; we are reporting the results obtained 
by discarding different numbers of low ^ values, obtained with the corresponding u value. 
The fits are generally of good quality, but the resulting continuum values of g^i depend 
noticeably on the number of skipped values. 

Our MC results are to be compared with the result of the form factor computation from 
Appendix D 

5r = 9.10(4). (6.9) 



40 



"'Skip 


i/K v-* ^ C! J 


rh i ^ /d of 




rhi^ /dnf 


1 
i 


O 1 0/1 ('fi'l 
y. iZ4(^0 ) 


o o /o 
z.z/ y 


O 1 00^'7^ 

y.izy(^ / j 


A A /Ci 

4.4 / y 


o 
Z 


o 1 onn 
i).iA\J[i\j ) 


O /I /s 
Z.4/ o 


O 1 O/l (Ci \ 

y.iZ4(^y J 


6.0/ o 


q 
O 




1/7 
A. i/ 1 


o 1 'io(^ o\ 

i).ioA[iA ) 


O.U/ ( 


/I 




1 /fi 


Q 

it . Lo lyio J 


O.U/ u 


5 


9.140(25) 


1.4/5 


9.148(22) 


1.6/5 


6 


9.134(33) 


1.3/4 


9.138(31) 


1.4/4 


7 


9.108(40) 


0.4/3 


9.112(38) 


0.3/3 


8 


9.099(43) 


0.8/2 


9.105(41) 


0.9/2 



Table 9: Fits of gR{K, oo) to a linear function of Q; nskip denotes the number of low ^ families discarded. 
The two columns correspond to the two rows of Table 8. 



6.3 The current correlation function 



Here we compare the bootstrap result for the current two-point function (2.11) with the 
lattice measurements. The extrapolation of the lattice data to the continuum is done 
by means of a two-step procedure, which is a variant of the method used for g^R. We 
first perform a FS analysis for those relatively short correlation lengths for which we 
could afford to measure the correlation function on lattices of large physical size. Using 
the FS scahng coefficients determined this way we are able to extrapolate the results 
of our measurements, corresponding to moderately large physical size, to infinite size. 
In the second step we take the continuum limit by extrapolating our results for infinite 
correlation length. 

For the first step we adopt an additive form of the FS scaling hypothesis: 



I{q;K,L)^I{q;K,oo) + z^'e 



C',{q;K) + Cl{q;K)- + 
z 



(6.10) 



We note that for the WARD case the subtracted correlation function is a linear combina- 
tion of a 2-point function and a 4-point function; therefore by the arguments used before 
one expects — A4 — 1/2. 

Compared to the case of g^,, the analysis here is complicated by momentum dependence 
and also the fact that we have used two alternative definitions (SUB and WARD) of 
the current correlation function. To be able to compare the results with each other and 
with the bootstrap calculation we used the dimensionless momentum variable q :— p/Mr, 
where Mr is the second moment mass. We interpolated our lattice results to integer q 
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values g = 1, 2, . . . , 50 by fitting tlie measured values with the 8-parameter formula 

10 2 

^(^) = E (6-11) 

fc=3 ^ 

This formula is motivated by the spectral representation and with the 8 parameters {ruk} 
it gives excellent representation of the current correlation function for all our lattices in 
the range q < 50. 

We fitted the leading coefficient CQ(g) (ignoring its K dependence and any sub leading 
terms) using the four data points of families 2, 4 and 6 for the SUB case and families 2 
and 4 for the WARD case. (Unfortunately the WARD data are not available for family 6, 
except for the lattice 6 with L = 200.) The results of these fits for q = 25 are shown in 
Figure 9. The fact that the linear fits are nearly parallel shows that our assumption of 
FS scaling works here. More importantly one sees that the difference between the two 
definitions disappears at large z, as it should''". We also note that FS corrections have 
opposite signs for the SUB and WARD cases and that the latter are much smaller. These 
qualitative features remain valid also for other g-values, although for small momenta the 
FS corrections for WARD data are no longer minute compared to the SUB ones at the 
same q. On the other hand for larger momenta the difference is even more pronounced. 

0.46 
. 45 
. 44 
0.43 
0.42 
0.41 
0.4 
0.39 

0.005 0.01 0.015 0.02 0.025 0.03 0.035 

Figure 9: Test of FS scaling: /(25) in SUB and WARD definition versus y^e"^. The three (approx- 
imately) parallel lines with negative slope correspond to the SUB data; increasing values correspond to 
families 2, 4 and 6 respectively. The two (approximately) parallel lines with positive slope are the WARD 
data for families 2 and 4. There are no WARD data available for family 6. 

Before turning to the infinite volume extrapolation let us briefly digress on the relative 
size of the statistical errors in the WARD and SUB data. In Figure 10 these errors are 

^Note that in this analysis we replaced ^ by 2;'. 
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shown as a function of the momentum q for both methods, for lattice 12. One sees that, 
with the exception of the first few points, the WARD data have much smaller statistical 
errors, for large momenta by about an order of magnitude! This remains true for all other 
lattices. The explanation is that before subtraction, the zero momentum component of 
the Fourier transform of the current correlation function has the largest fluctuation and 
its fluctuations decrease with momentum. Already at g ~ 4-5 the fluctuations are much 
smaller. If we subtract the zero momentum component, all the SUB points inherit its 
large error and beyond q > h the errors are practically constant. On the other hand, 
since the action density is known with a very good precision the errors of the WARD data 
points are almost the same as the unsubtracted ones and rapidly decrease with increasing 
momentum. For q less than 3 ~ 4, the SUB data have smaller errors since the fluctuations 
of I{q) and 7(0) cancel due to their strong correlation. Because of these two advantages 
of the WARD method (for not too low momenta), from now on we will use the WARD 
data exclusively. 




10 20 30 40 50 



Figure 10: Absolute statistical errors for SUB (increasing) and WARD (decreasing) current data as a 
function of q, for ^ = 421, L = 2100. The continuous lines are fitted for convenience. 

Let us now address the extrapolation to infinite lattice size. This is done by determining 
the FS coefficients in (6.10) from a small reference lattice and then use them to do the 
extrapolation for all other lattices. We used lattice family 4 to determine the coefficients, 
which is the family with four sizes and the largest coupling. Using family 2 instead leads 
to slightly different results that allow one to estimate the systematic error in the extrap- 
olation procedure. For illustration let us quote the (absolute) statistical error stat(g) and 
the (absolute) systematic error sys(g) obtained thereby at g = 5, 15, 25: stat(5) = 0.0002, 
sys(5) = 0.0006, stat(15) = 0.0001, sys(15) = 0.0003, stat(25) = 0.0001, sys(25) = 0.0001. 
One sees that the errors are small and the infinite volume extrapolation is under control. 

In the final step the extrapolation to infinite correlation length has to be performed. 
Since the largest lattice 12 already corresponds to a correlation length of ^ = 418 one 
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10 20 30 40 50 

Figure 11: Thermodynamic values of the current two-point function at ^ = 418 versus q. The solid line 
is the 2 + 4 particle form factor result and the dashed line is the absolute upper limit. 

might be tempted to regard this as superfluous. However if one were to take these data 
as representing the continuum limit one would have to conclude that the XY QFT does 
not coincide with the 0(2) bootstrap theory! This is because, in contrast to the 
measurement, the statistical errors here are very small and the data differ significantly 
from the bootstrap result. Moreover, as explained below, the truncation error in the 
bootstrap computation is under good control. The situation is illustrated in Figure 11. 

In general it is difficult to strictly control the systematic error in a form factor computation 
caused by the truncation in the number of intermediate particles. One only knows that the 
truncated result provides a strict lower bound on the exact answer, since (for a two-point 
function) all multi-particle contributions are positive. In the 0(2) model however we also 
have a strict upper bound. This is because the exact I{q) is known to be increasing and 
to approach the value 2/7r ~ 0.637 at infinite momentum; c.f. (2.17). On the other hand 
the 2+4 particle approximation is likewise increasing and approaches 0.621 at infinity; the 
difference 0.016 is an upper bound on the error made by the truncation, because also the 
higher particle contributions are monotonically increasing. The true value of the bootstrap 
function I{q) is somewhere between the 2+4 approximation and this approximation + 
0.016, shown as a dashed line in Fig. 11. For the relatively low momentum range we are 
interested in, the true value is probably closer to the 2+4 value than to the upper limit. 

If we use all our data to perform an extrapolation to infinite correlation length, the 
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situation changes drastically. We assume a cutoff dependence asymptotically linear in Q 

I{q; K, oo) = 7(g) + I'{q) Q + 0{Q') , (6.12) 



and extrapolate our measurements to infinite correlation length by means of a linear fit 
to the data points corresponding to families 10, 11 and 12. Our thermodynamic data 
(for families 2,4,6,9,10,11,12) together with this fit is shown in Figures 12 and 13 for 
g = 5, 15, 25. For concreteness we used u = 1.33, corresponding to ngkip = 4 in Table 
6, in all our fits. Varying u around this value we observed that our results are rather 
insensitive to the precise choice of u in the range 1.2 - 1.5. 
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Figure 12: Thermodynamic values of the current correlation fimction for lattice families 2,4,6,9,10,11,12 
for g = 5 versus Q. The solid horizontal line shows the 2 + 4 form factor result and the dashed horizontal 
line is the absolute upper limit. The linear fit is based on the three biggest lattices, corresponding to 
families 10, 11 and 12. 

.41 ^- - _ . _ 

0.405- 

0.40- 
0.395- 

0.39- 
0.385- 

0.38- 

. 375 L1 .... I .... I .... I .... I .... I .. . 

0.05 0.1 0.15 0.2 0.25 




0.46- 
0.45- 
0.44- 
0.43- 




. . . I .yj . 42U1 , , , , I , , , , I , , , , I , , , , I , , , , I , , , , I , , , , I ,• 

.3 0.35 0.05 . 1 0. 15 .2 .25 .3 0.35 



Figure 13: Left: same as Fig. 12 for q = 15. Right: same as Fig. 12 for q = 25. 

From these figures one infers that the approach to the continuum limit, for g > 10, 
is non-monotonic. While this is not unusual in itself, it is most remarkable that for the 
momentum range g > 20, the turning point is around C, ~ 40, a rather large value. Beyond 
the turning point the data behave as expected and follow a curve that is approximately 
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linear in Q. If wc extrapolate our measurements to the continuum limit using the linear 
fit, the extrapolated points agree reasonably well with the form factor calculation. This 
is shown in Figure 14, which is one of the main results of this paper. According to this 
figure the current two point function of the XY QFT is very close to that of the 0(2) 
bootstrap theory. The absolute difference between the extrapolated points and the 2+4 
FF result is in the range [0.0003,0.0032]. This difference is positive and much less than 
0.016, consistent with the hypothesis that it is due to higher particle contributions. 
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Figure 14: Thermodynamic values of the current two-point function extrapolated to the continuum 
limit versus q. The solid line is the 2+4 particle form factor result and the dashed line is the absolute 
upper limit. 

Not only the extrapolated continuum values I{q) but also the coefficient I'{q) governing 
the rate of approach can be compared to theory. In the theoretical framework presented in 
ref. [4] (and recalled in Section 4) this coefficient can be calculated from the continuum SG 
theory. In [4] this was done in asymptotically free perturbation theory at two-loop order. 
This is expected to be vahd for large momenta subject to the additional constraint [4] 
\ogq < v^3/2g. Taking Q ~ 0.1 our data for 10 < g < 50 are just in the window 
of validity. Here we calculated the leading 2-particle contribution to the first correction 
coefficients using the known [24] 2-particle form factor of the Noether current in the SG 
theory. In the lattice theory we define I'{q) as the slope of the linear fit on the data 
for the 3 largest lattice famihes, 10,11,12. (Recall that all our results correspond to the 
choice u — 1.33.) The comparison with the theoretical results is shown in Figure 15. The 
agreement is quite remarkable (for the value of u chosen). In particular both analytical 
computations predict a change of sign in I'{q) between q — lQ and q — 20, which is indeed 
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observed for the fitted numbers obtained from the MC measurements. 






2 















1 


- 




/ *2 


, 1 , , , , 1 














30 40 


50 





1 















2 















3 













Figure 15: Fitted values of the coefBcients /'(g) versus q. The sohd curve is the perturbative result of 
[4] and the dashed curve is the leading (2-particle) form factor contribution. 



6.4 The spin correlation function 

Finally we consider the two-point correlation function of the spin operator in Fourier 
space. While the current correlator had to be subtracted, the spin correlator has to be 
normalized (multiplicatively) before it can be compared to analytical calculations. We 
define 

G<-.(,) ^ |M. (6.13) 

where G^i^q) is the bare spin correlator. Traditionally one takes /i = 0, which amounts to 
normalizing the spin correlator with the susceptibility x- A problematic feature of this def- 
inition is that, just like with the current defined by the SUB method, the statistical errors 
are big, dominated by the large error of X. We thus introduced G^^^{q) := {qj \xfG^^\q), 
whose values are closer to unity, and decided to take /j, — 5. For q > 5 this leads to much 
smaller errors, in many cases smaller by a factor 3-4. We adopted this unusual choice 
of normalization because in this paper we focus our attention to the range q > 5. The 
reason for concentrating on this range is that it is here the cutoff effects show interesting 
non-monotonic behaviour. (For the low momentum range ^ < 5 it would be better to 
take the normalization at /i — 0, but this is not investigated here.) 

Adopting the normalization at g = 5 the FS analysis is analogous to the current case. 
The FS scahng hypothesis is apphcable and allows one to extrapolate G^^\q) to thermo- 
dynamic lattices. The results are shown in Figures 16, 17 and 18 for q — 15, 25 and 35, 
respectively. Again, similarly to the current case, one sees a non-monotonic approach to 
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the continuum limit, with a turning point which is for g > 20 around ^ ~ 100. The actual 
points are still significantly away from the analytical prediction, but beyond the turning 
point they move into the right direction. We did not attempt to fit a linear function to 
the data but our results for the spin correlator are not inconsistent with the form factor 
bootstrap. 



1.26 
1.24 
1.22 
1.2 
1 . 18 
1.16 
1 . 14 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 

Figure 16: Thermodynamic values of the spin correlation function &^\l5) for lattice families 2, 4, 6, 
9, 10, 11 and 12 versus Q. The dashed line shows the 1 + 3 form factor result. 
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Figure 17: Thermodynamic values of the spin correlation function &^\25) for lattice families 4, 6, 9, 
10, 11 and 12 versus Q. The dashed line shows the 1 + 3 form factor result. 
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Figure 18: Thermodynamic values of the spin correlation function &^\35) for lattice families 6, 9, 10, 
11 and 12 versus Q. The dashed line shows the 1 + 3 form factor result. 

7. Conclusion 



Since we already surveyed our motivation and some of the theoretical issues involved in 
the introduction, let us return here to the question raised in the title. Screening the 
comparison between bootstrap and lattice theory for the various quantities considered, 
we would tend to answer the question in the affirmative. Probably the strongest Pro 
argument stems from the intrinsic coupling g^. The final results in both approaches have 
an estimated (systematic) error of less than one percent, so that the good agreement is 
remarkable. For the two-point function of the Noether-current the prediction [4] for the 
lattice artifacts could be tested. After, and only after, the lattice artifacts are taken into 
account a good and non-trivial agreement with the form factor result emerges. A final 
decision whether the remaining small differences are due to the neglected higher particle 
contributions or signify in fact a true difference of the two constructions cannot be reached 
at this stage. Quantitatively the least convincing are the phase shift results. However in 
the lattice framework measurements of the phase shifts are technically difficult and here, 
as in other models, mainly the qualitative features at low energies can be probed. But 
the latter do agree with the bootstrap prediction. Each comparison considered separately 
certainly leaves room for doubt, but collectively they do suggest that the continuum limit 
of the XY model and the 0(2) bootstrap theory are the same QFT. 

Concerning future work, a more detailed exploration of the superselection structure should 
be interesting. The new parafermionic superselection sector found here is probably accom- 
panied by a third ('disorder-like') sector. Their interplay e.g. on the level of the operator 
product expansion as well as an explicit field theoretical construction remains to be found. 
Finally, as a test case for other sigma-models, it would be important to understand which 
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quantities in the XY QFT can be understood, qualitatively or quantitatively, in terms of 
a perturbed conformal field theory description. 
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A. Quantum group covariant form factors 

Here we derive the necessary and sufficient conditions (3.13) on the statistics phases that 
ensure the quantum group covariance of the form factors. The additional conditions hke 
(3.27) required for multiplets transforming irreducibly are also detailed. 

To fix our conventions we begin by recalling a few definitions for the action oiUq{su{2)) on 
some irreducible representation. Background material on quantum groups in 2-dimensional 
physics can be found in the book [18]. The Hopf algebra Uq{su{2)) is generated by X±, H 
that act for generic q on an irreducible spin j module according to 

X±\j, m) = ^J [j T m]q\j ± m + l]q \j, mil), 

H\j, m) — 2m |j, m) . (A.l) 

Here |j, m), m G {— j, — j + 1, . . . , j — 1, j} denotes a basis of the {2j + 1) -dimensional 
irreducible module 2j + 1, and [n]q = (g" — q^'^)/{q — q^^)- For q — — e~*'^/^, p > 3, an 
upper bound on the allowed isospins j exists. It reads j < p/2 — 1 and is related to an 
enlarged center; see e.g. [2]. To the best of our knowledge the case q — —1 has not been 
studied explicitly in the literature, but it is not hard to work out the aspects needed here. 
First, as a Lie algebra U-i{su{2)) is isomorphic to su{2), but the co-multiphcation differs 
by signs. Guided by the sample computations presented below and the formal p — > oo 
limit of the above relation, we expect that for = —1 no trunction of the allowed isospins 
occurs. For definiteness we fix the roots q^^"^ — i, q"^!'^ — —%. 

For generic q the comultiphcation of Uq{su{2)) is 

AH = H + , (A.2) 

for g — s> — 1 we define it with the above choice of roots. For j = 1/2 we write |i) = 
|l/2, il/2) for the basis of the defining representation 2. In the n-fold tensor product 
2*^"" wc write |a„,...,ai) := |q;„) Cg) . . . ® |ai), o^j G {i}) ior the natural basis. The 
'charged' components of a form factor are introduced as the coefficients with respect to 
this basis, i.e. 

I/) = XI f»r....ai\an---ai)- (A.3) 

«nv,«l 

By construction the quantum group generator H'^"' is diagonal on this basis and its 
eigenvalues e := «„ + ... + ai are the U(l) charges used in section 2. The raising and 
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lowering operators A^"'>X± act as a 2" x 2" matrix S-t on the basis \an, ■ ■ ■ ,Oii). We 
choose a lexicographical ordering of the basis vectors that is symmetric under the flip 
aj — > — «j. Then S_ = S^. Further there is an induced action on the coefficients fa„...ai 
in (A.3) which is implemented by S_ for A(")X+ and by S+ for A^'^^X.. The 2" x 2" 
matrices E± are triangular and 'sparse' with only a few blocks different from zero. The 
block structure arises because evidently A(")X± maps the charge e sector into the charge 
e ± 2 sector. Explicitly the matrix E_ acts on the form factor components as 

S- : /a„...a. e^^(-^) ^ (-1)"-Va„........... , (A.4) 

and similarly for E+. The mapping (3.14) could be used to 'untwist' the 5'C/_i(2) co- 
multiplication, i.e. to remove the phases in (A.4). We refrain from doing so because the 
'untwisting' does not induce a physically interesting correspondence between the form 
factor sequences of the SU{2) and the SU-i{2) bootstrap theories. 

As usual the n-particle form factors carry a representation of the permutation group Sn- 
Its representation matrices are 2" x 2" matrices Ls{9), s e Sn- The quantum group 
invariance of the S-matrix (3.7) generalizes to 

E± L,{e) = L,(^) E± , ^seSn. (A.5) 

For completeness let us note the explicit definition. One sets 

LsM7.:t = €:■■■ s^f^^iiej^i,) . . . C > j = i, . . .n - 1 , (a.6) 

for the generators Si, . . . , s„-i, acting by Sj{9n, . . . ,9i) = {On, .... 6j, ^^j+i, . . . , 6*1) on the 
rapidities 9 := {9n, ■ ■ ■ ,0i). By means of Lss'id) = Ls{9)Lsi{s^^6) this extends to all 
s, s' e Sn- The invariance (A.5) clearly entails that the form factor Eq. (3.9a) (and its 
generalization to generic s G S^} are covariant under the quantum group action. 

It is natural to ask whether the same can be achieved for the cyclic form factor equations 
(3.9b). In that case the cyclic Eqs in the charge e sector 

fa„...aiiOn + 27li, 9n-l, ■ ■ ■ , Oi) = r]a„{e) /a„_i...Q;iQ;„ (^n-1, ■ • • , 6'„) , (A.7) 

with e = a,„ + . . . + «! and those in the charge e ± 2 sector would again be compatible 
with the quantum group symmetry. Explicitly this means that starting from (A.7) and 
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performing e.g. the substitutions (A. 4) the result should be an identity by virtue of the 
cyclic Eq. in the charge e ± 2 sectors. This condition gives rise to a set of overdetermined 
relations among the phases Tia{e), - which turn out to be self-consistent. Using (A. 4) the 
solution (3.13) can be verified. 

Finally we turn to the residue equations. Consistency requires that the inverse of the 
matrix in (3.9b) appears on the right hand side, irrespective of its concrete form. In 
the charged basis and for n > 3 one has: 

-TCS0^,=e„_i+iTTfa„...ai {On, ■ ■ ■ , Ol) (A.8) 
= k(e)-^^.„-....s,(^)£:f«:-_t-. - ■ ■ ■ C] U.-...M0n-2, . . . , ^l) . 

Here e — an + ■ ■ ■ + ai — Pn-2 + ■ ■ . + Pi refers to the charge sector. For the specific choice 
of phases (3.13) these Eqs can be seen to be likewise quantum group covariant. 

In summary there exists a preferred (and up a normalization uniquely determined) choice 
of the statistics phases i]a{e) for which the form factor equations (3.9), (A. 9) are covariant 
with respect to the quantum group W_i(sm(2)). This means its solutions can be grouped 
into multiplets that transform covariantly under the symmetry group, and one can restrict 
attention to those transforming irreducibly. Technically the irreducibility condition can 
be encoded into a parameterization of the form factors that is adapted to the embedding 
of the irreducible spin j module 2j + 1 into 2®", including multiplicities. Essentially it 
amounts to determining the generalized Clebsch-Gordon coefficients. 

To faciliate the comparison with the familiar l{i{su{2)) = su{2) case we first consider 
the decomposition for generic q and specialize to g = — 1 only at the end. Thus let 2 
again denote the defining two-dimensional representation of Uq{su{2)) and consider the 
decomposition of 2*^" into irreducible representations. It assumes the familiar form 

2««= m,(n)(2j + l), (A.9) 

jo<j<n/2 

where mj{n) is the multiplicity with which 2j -|- 1 occurs and jo = 0, 1/2 for n even,odd, 
respectively. For generic q these multiphcities are the same as for su{2), only the Clebsch 
Gordon coefficients differ. As outhned before we expect the hmit — > — 1 to be regular in 
the sense that no trunctation of the isospins occurs and that the multiplicities mj{n) are 
the same as in the undeformed case. The multiplicities then are conveniently computed 
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from a generalized 'Pascal triangle' described by the recursion relations 
mn/2{n) = 1 , mj{n) = , j <0 , 

mj{n) = mj_i/2(n - 1) + mj+i/2(n - 1) , j = jo, . . . ,n/2 . (A. 10) 

The highest weight conditions A^^'^X^v = 0, v E 2®", arc readily solved and yield Tnj{n) 
linearly independent solutions on each of which a spin j multiplet can be based. 

Below we hst for n = 2,3, 4 a basis for the spin j sector in 2®". The multiphcites are 
taken into account by displaying a mj(n)-dimensional family of highest weight vectors. 

n = 2 : 

i = 0: q'''\ + -)-q-'''\-+), 
n = 3 : 

i = l/2: ^;l/2 = All -++) +A2I + -+) +A3I ++-) , g^Ai + gA2 + A3 = , 

= (A2 + g"^Ai)| + — ) + (Ai + As)! - +-) + (5A1 + A2)| - -+) , 
j = 3/2: E^|+++), A; = 0,1, 2, 3. 

n = 4 : 

i = : vo^\-- + +)-{q- q-')\ - + -+)- | - + + -)- | + --+) + + + - -) , 

Vo = q~'\- + -+) + \- + + -) + \ + --+)-q\ + - + -) , 
j = 1 : = All - + + +) + A2I + - + +) + Aal + + -+) + A4I + + + -) , 

g^Ai + g^A2 + ^Ag + A4 = , 

E_vi ~ g(gAi + A2)| - - + +) + g(Ai + Ag)! - + -+) + (A3 + gA4)| - + + -) 
+ (gA2 + Aa)! + --+) + (A2 + A4)| + - + -) + (Ai + 5-^4)1 + + --). 

E'_v, = g^Ail + ) + q'[X2 + {q- q-')Xi]\ - + - -) 

+ [A3 -(g_5-i)A4]| -- + -)+ A4I +), 

i = 2: Ei|+ + + +), A; = 0,l,2,3,4. 

To illustrate the use of this table let us look at the n = 3, j = 1/2 entry. The form factor 
components /i, /2, /s in (3.23) play the role of the A's and for q = —1 one obtains (3.27). 
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B. 3-particle Form Factors 



Here we construct the 3-particle form factors of the spin and the parafermion field by an 
adaptation of the technique of [3]. From the SG viewpoint these fields are nonlocal which 
is why they have not been considered in [3]. All even particle form factors of the SG fields 
were constructed by Smirnov [43] where also the Bethe ansatz technique instrumental in 
[3] is implicit; see also [47] for related results in the mathematical literature. 

Adapting theorem 4.1 of [3] the 3-particle FF can be represented as the contour integral 

JC (B.l) 

Here 

1^(^3, ^2, ^i) = y{e^ - 62) y{es - e,) y{92 - 61) , (B.2) 
and the only nonvanishing component of Tq^q.^ is 

where 

^ ^-m^-i-s .^y(- + —\t(-—]. (B.4) 

Finally the "pseudo- vacuum" vector is 

The integration contour C consists of several pieces. First, there are three small clockwise 
circles around the three points ^1,^2 and ^3. In addition C also contains a line integral 
parallel to the real axis such that the integration path goes between 6m — ^tt and 9m — 2i7r 
for all the three 9m- 

Since C is defined relative to the arguments 9m, when we analytically continue (B.l) it is 
best to deform the contour parallel to the arguments. This way it is trivial to see that 
the solution satisfies (3.8). It is also relatively easy to see that the Bethe Ansatz hke 
construction (B.l) ensures that (3.9a) is also satisfied, independently of the contour C. 
(To show this one has to use the Yang-Baxter equation satisfied by the S-matrix.) It is 
more difficult to verify (3.9b) because here one has to take into account that the contours 
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are different on the two sides of tfie equation. Similarly for the residue equations (3.25) 
the contour is different from the original C. 

Inserting the S-matrix (3.1) Eq. (B.l) can be rewritten as 

" 3 



fm{03, 02, = AfY{es, 02, e,) du e^(<^i+^2+(^3-2«) 



tm{03,92,ei;u), (B.6) 



where 



t2{9s,92,9i;u) = 
h{9^,92,9i;u) = 



9?.-u 9o-u 



in 



iTT — 9s + u iTT — 92 + u iir — 9i + u ^ 

9s — u m 

m — 9s + u m — 92 + u'' 

in 



(B.7) 



in — 9s -\- u 



Next we examine whether the solution (B.7) is compatible with SU_i(2) symmetry, which 
requires the vanishing of the linear combination (3.27). This can be written as 



C(^3, 92, 9,) = ^fY{9s, 92, 9,) [z,{9s, 92, 9,) + ^2(^3, ^2, ^1)] , 



(B.8) 



where 



zi{9s, 92, 9i 



Z2{9s,92,9i) 



Using the identity 



(B.9) can be rewritten as 



+e2+dz-2u) 



9k -u 
in — 9k + u 

k=l 



n 

.A;=l 
3 



(f){9k - u) 



ITT — Z 



(f){z) — 0(2; — 2m) 



(B.9) 
(B.IO) 

(B.ll) 



^1(^3,^2,^1) = e^-- / , 
Jc+ 



^s(0i+6>2+e3-2«) 



.k=l 



(B.12) 



where the contour C+ is shifted by 27ri, i.e. it consists of three small circles around 
9m + 27ri and the line integral goes between 9„i + ^tt and 9.„i. From (B.4) we can see that 
the small circles do not contribute here since the integrand is regular there. The only 
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relevant singularities are those at 6m and it is easy to see that the contribution of the 
shifted line integral is precisely the same as that of C. We thus have 

C(^3, 02, e,) = + i)MY{es, 62, e,)z2{es, 62, e,) . (B.13) 

This proves the assertion after (3.27). Remarkably the quantum group invariance is not 
visible on the level of the integrand in (B.l) but only after the integral has been performed. 
In addition this fixes the value of the spin to be s = ±1/4, without mod(l/2) ambiguities. 
This is because the integral in (B.l) exists for \s\ < 3/4 only. 



C. 4-particle Form Factors of the Noether Current 

The form factors of the Noether current can be found in Smirnov's book [43]. We have 
adapted this result to our notations and conventions for the 4-particle case. 

Let us introduce the reduced form factors g that are defined by 

(01,92,9^,94) ^Y{0,,02,03,04)9 [01,02,0^,0^), (C.l) 

where 

Y{0,,02,0M = \{y{Oi-0j). (C.2) 

Using the 0(2) symmetry and charge conjugation, we need only the following components: 

9++ — (^1,^2,^3,^4) = —g — ++(^1,^2,^3,^4), 

9+ — +(^1,^2,^3,^4) = — 5'-++-(^i, ^2, ^3, ^4) , (C-3) 

g+-+-{9i,92,9z,9i) — —g-+—^.{9i, 92, 93,94) . 

Further, using the axioms, we can express everything in terms of a single function 
M^i, 92, 93, 94) as follows. 

g++ — (^1,^2,^3,^4) = ^(^1,^2,^3,^4), 

g+^^49,,92,93,94) = ^(^4 + 271^,^1, ^2, ^3) (C.4) 
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and 



^2, ^3, ^4) = U{9,, 03 + 2711, 9,, 62) 

"4 — "3 I 

?7r "I 
- -. ^—w- M^i, ^4 + 27ri, 0s, 62) . 

ITT — O4 + O3 J 

This function is given by 

2^e-2A(o) -i(E,ti/3i) 

where 

2 4 

p{(3,,(32,(3s,(3,) = (A + /32 -(3s- (34- 3. _ 3. _ 

1=1 j=3 * 

and 



/(/?i,/?2,/?3,/34) = / dae°'{q{a,pi)q{a,p2)x{a,ps)x{a,p4) 
J —00 

+ a;(Q;,/3i)a;(Q;,/32)g(a,/33 + 2m)q{a,(3A + 27ri)} 

Here we defined 
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D. Calculation of the subleading correction to 

The leading contribution to the intrinsic couphng in the 0(2) model was calculated in 
[7]. Here we present the calculation of the most important subleading term reducing the 
theoretical uncertainty in this quantity to a few per mille. This appendix relies heavily 
on [7], especially Sections 3, 6 and Appendix C. We use the notations and conventions of 
that paper. 

The first few contributions to 74 for the Ising model and for the 0(3) model were also 
calculated in [7] . Let us compare the results. 



contribution 


Ising model 


0(3) model 


121 


-4.99343 


-4.16835 


123/1 


-0.01348 


-0.01351 


123/2 


0.10610 


0.11901 


123/3 


0.00000 


-0.00200 


141 


-0.00265 


-0.00407 



Table 10: The first few contributions to 74. 123/i stand for the contribution of the integrals V^'^^ for 
i=l, 2 and 3. 



The pattern is strikingly similar for the two models. For the XY model we so far only 
have the 121 contribution, which is somewhere inbetween the corresponding values for the 
Ising and the 0(3) model. If we assume that the corrections follow the same pattern also 
for the 0(2) model, already the calculation of the 123/2 term yields 74 with a precision 
better than one percent. We will see that we have all the ingredients necessary for this 
calculation. 

Using the FF axioms and some of the expressions in Appendix C of [7] we have 

g^^\(3,ai,a2) = gl^,,{i7r, (3, -(3) Gl^^^^iin - (3,ai,a2) gtl^^^{(3,ai,a2) (D.l) 
- ^i66(«7r, -/3, (3) GL^^^im - P, ai, G*l^^bi(^i, 0i2, (3) . 



Now it is easy to see that g'^^-*(0, a, —a) = so 



V^(^) = 



(D.2) 
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in general. Further using (D.l) in {CAS) of [7] we get 

- GUtn, -/?, /?) [gtU-P, 0, /?) + giliO, -/?, /?) } . (D.3) 



We see that 



so also for generic n 



y(^) = . (D.5) 
(D.2) and (D.5) together imply that the conjecture (C57) of [7] is true. 
Similarly simplifying (C.37) of [7] gives 

9^'\P) = GU^^, P, -P) OmiP, -P, 0) (D.6) 
and using this in (C.36) of [7] we have 

^i\,(i7r, V, -v) J^U-V: V, 0) . (D.7) 
It is easy to see that we already have everything that is necessary to compute (D.7) since 

Kyli-V, V, 0) = Syi,gp{v) f^g^iiTT + v) . (D.8) 

Putting everything together we have 

V^'^ = ^ Tdv ^ e^(^) M{v) , (D.9) 
D47r Jo cosh V 

where 

H{v) = A(i7r + v) + A(i7r - v) + A{2v) + A{-2v) + A{v) + A{-v) - 2A(0) (D.IO) 
and 

M{v) = -r-^ — ^ ^\K{v)K{i7r + v)\27r'^ + n{n-2)v{i7r -v)] 

[in — V) [{n — 2)v — 27ri\ 

+ 2K{v)L{in + v) [{n + + {n - 2)v{iTT - v)] (D.ll) 
+ 2L{v)K{i'K + v) [27r^ + (n - 2)v{ iir — v)] 
+ AL{v)L{i7T + v) [27r^ - iirv + (2 - n)v^] j . 
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To summarize, using the results of [7] (D.9) in the n = 2 case can be represented as 



36 0077^ Jo 



d9 



sinh^ 9 
cosh*^ 9 



e^+^3+i3i ^(4^2 + 257r2)(4^2 + 49^2) 



+ e 



{49' + 97r^) 



^ 2(e2 + 7r2)J 
(D.12) 



where 



H(9) 2 / coso;^ + cosh7ra;(cos2a;^ + coso;^ — 1) — 2 
^ ^~ I CO sinh7ra;(l + e'^'^) 



Dm 

Ds{9) 
Dm 



duj cos - 1 
V 2 sinh ^ 



A;(a;) , 



do; cosh TTo; cos ^ — cosh ^ , ^ 



sinh TTo; 



du; cos 'f - cosh , , , 
UJ smh TTo; 



Further 



and finally 



where 



+ iA2 = i{i7r - 29){3i7r + 29){5tn + 29)e'^^^^^ 



Dm = - — 



^ dcu sin 



2 



Numerically we find 



UJ 2 cosh ^ 



y(2) ^ 0.00724518. 



k{uj). 



(D.lSa) 
(D.lSb) 
(D.13c) 
(D.lSd) 

(D.14) 

(D.15) 

(D.16) 
(D.17) 



Now we are in a position to be able to fill in some 0(2) entries in the previous table. 

As expected, the available 0(2) data follow the same pattern as before. Furthermore the 
0(2) numbers are in between the 0(1) and 0(3) ones. So (with some confidence) we can 
predict the uncalculated (*) contributions to be close to the avarage of the corresponding 
Ising and 0(3) entries. In this way we get 



(*) = -0.01786 ±0.00893, 



(D.18) 
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contribution 


Ising model 


moQei 


(J[o) mociei 


1 91 


— ^t.yyo^to 


— 41.00 ( io 


— 41. iOOOO 


123/1 


-0.01348 


* 


-0.01351 


123/2 


0.10610 


0.11592 


0.11901 


123/3 


0.00000 


* 


-0.00200 


141 


-0.00265 


* 


-0.00407 



Table 11: The first few contributions to 74. 123/i stand for the contribution of the integrals F^'^ for 
i=l, 2 and 3. 

where (generously) we allowed for 50% error here. This gives a total k + I + m — 6 
contribution of 

0.09806 ±0.00893. (D.19) 

Estimating the k + I + m > 8 contributions, as usual, to be less than 10% of (D.19) our 
final estimate is 

74 = -4.559 ±0.019. (D.20) 

We also need the product 7252- We have computed numerically the 3-particle contribution 
to both 72 and S2 using the 3-particle form factors constructed in Appendix B. We found 

72;3 = 0.001813(1) and 62-3 = 0.00003016(2) . (D.21) 

Thus 72^2 = 1.00184 and finally we get with an error of a few permille: 

5r = 9.10(4). (D.22) 
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E. Test of random number generators 

Since wc are interested in achieving numerical precision for many quantities to an accuracy 
of < 1%, a considerable source of concern to us was the random number generator (RNG). 
Indeed at an initial stage of this project we found large standard deviations between results 
obtained by various generators. 

Our first test concerned comparison of computations using various RNG's, with exact 
results on small lattices. The (practically) exact result on a 3 x 3 lattice is obtained by 
discretizing the spins, taking 0(2) — > Z(A^) and summing over all A^^~^ terms* in the 
partition function. The convergence to the 0(2) case is extremely (exponentially) fast. 
As illustration in Table 12 we give the values of the susceptibility for K — 0.25, L — 3 
and N — 6, . . . , 10. Some generators already failed this test. The exact numbers were 
also useful to check our programs. 



N 


X 


6 


1.7619848372 


7 


1.7619804581 


8 


1.7619803546 


9 


1.7619803525 


10 


1.7619803524 



Table 12: The susceptibility of the Z(A^) model on a 3 x 3 lattice at K = 0.25. 

As a next step we compared results obtained by different RNG's on larger lattices, see e.g. 
Table 13 where we tabulated our results for the susceptibility at K — 1.0, L — 256. Here 
rand is the RNG hsted in Language Reference XL Fortran for AIX (Version 3 Release 2) 
and SGI is the RNG provided by Silicon Graphics for the SGI 2000 machine. The nag (the 
g05caf RNG by Numerical Algorithms Group) and ranlux [34] are portable RNG's. The 
latter has a single- and a double precision version, both with an extra choice, a "luxury 
parameter" . The notations "rbcs-O" , "rbcd_l" and "rbcd_2" refer to the precision and the 
value of the luxury level parameter. 

To our knowledge ranlux is the only generator known with proven randomness qualities. 
Unfortunately, for "historical reasons" we used it only in the later stages of the project, 
while most of other runs were using nag. Reassuringly we found in all our tests, that the 

*Due to the global symmetry one spin can be fixed. 
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program 


KiNLr 


X 


i 


rand 


iouy.o^y ) 


1 


SGI 


1607.3(9) 


1 


nag 


1604.7(8) 


2 


nag 


1605.5(9) 


2 


rlxd_2 


1604.1(9) 



Table 13: The susceptibility at X = 1.0 and L — 256 using different RNG's and two 
different programs 

nag generator produced results consistent with ranlux. The combined nag resuh (same 
RNG but different programs) in Table 13 is 1605.1(6), which is only 1-sigma away from 
rlxd_2. Note however, the rand result is 4.2-sigma away from rlxd_2 while the SGI result 
is 2.5-sigma away. 

Although the latter deviation is still not too serious, the SGI RNG gave also suspicious 
results on lattices with very large physical size, z ~ 14: the data obtained by this RNG 
were several sigma too high above the FS scaling lines. Since the z > 10 data were not 
needed in the extrapolations anyhow, we simply omitted these data points from Table 4 
and did not use them in our fits. Nevertheless, to make sure that the discrepancy was 
indeed due to the failure of the RNG we remeasured a few of these points with ranlux. 
Table 14 shows the results for these two RNG's, together with the fits of Table 5. Note 
that the SGI results show large deviations from both ranlux and the FS fit, especially for 
the susceptibihty, which is always too high, by 5.6 to 7.6 standard deviations. 

The only SGI data present in Table 4 arc the 4 points aX K = 1.04 and 1.05 for z = 4 
and 5. These data, in contrast to the z = lA points, agree with the FS fits. We have also 
rechecked the K = 1.04, L = 578 point with rlxs_0 and got ^ = 142.10(9), x = 14142(12), 
gji = 7.604(9), in good agreement with the SGI results^ in Table 4 therefore we did not 
repeat all these measurements with ranlux. 

In one of the programs we measured the quantities with the standard estimator along with 
the improved one, and checked that they agree within the errors. The other program used 
a Ward identity for checking. Note, however, that the agreement in these quantities does 
not guarantee yet the correctness of the results: the error of the standard estimator is 

^Note, however, that at these very large correlation length our relative errors are much larger than 
those at smaller ^. 
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K 


r, 

1j 


c 
s 


■A/ 

A 


i/R 




0.92 


150 


1A r'r'A/o\ 

10.559(3) 


163.110(36) 


8.879(25) 


rand 






iA r'r'Ao/A\ 

10.5508(9) 


162.843(13) 


8.864(8) 


nag 






10.5499(7) 


162.829(10) 


8.869(6) 


rlxs_0 






10.5507(9) 


162.835(12) 


8.871(8) 


rlxd_l 






lU.t)OlU(ll ) 


1 o /I n /'I p.\ 
loz.o4U(lo ) 


Q 001 /'I r\\ 
o.ool(lU ) 


-.^1 O 

rlxa_z 






10.549(1 ) 


lo2.oo(2o ) 


O.O M (o) 


IPC fi + 

r o nt 


0.97 


OA A 

dOO 


21.o59(o) 


557.41(15) 


8.958(25) 


SGI 






O 1 POT/ A \ 

21.o27(4) 


556.32(11) 


8.952(6) 


1 o nt 


1.0 


5o0 


40.1oo(lo) 


1615.68(48) 


A AOT/ 0A\ 

9.037(29) 


SGI 






/I n 1 nv/ 1 o\ 
4U. 10/(12) 


1/^11 / /I o\ 

loll. /8(4o ) 


9.010(2o ) 


rlxa_2 






40.096(5) 


1611.4(3) 


8.989(6) 


FS fits 


1.02 


1000 


69.647(23) 


4184.5(1.2) 


9.066(34) 


SGI 






69.533(24) 


4174.5(1.3) 


8.962(33) 


rlxd_2 






69.505(20) 


4173.0(1.4) 


9.026(8) 


FS fits 



Table 14: The correlation length and susceptibility for z ^ 14 for different RNG's and 
the corresponding value obtained by finite size (FS) extrapolation taken from Table 5. 

usually much larger than that of the improved estimator, while even the bad z — 14 
results passed the WI test. 
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F. Vohwinkel's results for the 2-particle energies 

We here reproduce the original data table of Vohwinkel, giving the 2-particle energies 
obtained nearly 10 years ago [49]. We do not know the reason why Vohwinkel did not 
publish his results, but it could be that he just confused the assignment of quantum 
numbers, and so could not match his results with the proposed S-matrix. With the 
correct identification his data are listed in Table 15. 

As mentioned in Section 5 his values for the single particle masses on the lattices with 
z >~ 10 are in good agreement with ours. It is only on the lattice with K — 0.97 and 
L — 128 that our values differ by ~ 4 standard deviations. 



K 


0.97 


0.97 


0.92 


0.86 


L 


128 


256 


128 


64 


m 


0.04633(2) 


0.04620(1) 


0.09465(3) 


0.1711(1) 


2 


0.1019(2) 


0.0956(2) 


0.1949(3) 


0.353(1) 




0.1612(3) 


0.1149(3) 


0.2334(4) 


0.435(2) 






0.1481(6) 


0.2982(7) 


0.566(3) 






0.1872(8) 


0.3760(10) 




1 


0.1298(1) 


0.1038(4) 


0.2116(4) 


0.389(1) 






0.1324(4) 


0.2676(7) 


0.510(2) 






0.1705(6) 


0.3440(20) 







0.1081(15) 


0.0967(3) 


0.1966(10) 








0.1220(8) 


0.2431(20) 








0.1614(14) 


0.3143(40) 





Table 15: Masses and energies obtained by Glaus Vohwinkel in 1992 
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